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PREFACE 


This  report,  prepared  by  Gregory  B.  Baecher  of  NEXUS  Associates,  Wayland, 
Massachusetts,  with  assistance  from  D.  DeGroot,  C.  Erikson,  and  A.  Pais,  under 
Contract  No.  DACW39-83-M-0067 ,  provides  details  for  the  statistical  analysis  of 
geotechnical  engineering  aspects  of  new  dam  projects.  It  was  part  of  work  done 
by  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  in  the  US  Army  Civil 
Works  Investigation  Study  sponsored  by  the  Office,  Chief  of  Engineers,  US  Army. 
This  study  was  conducted  during  the  period  October  1983  to  September  1985  under 
CWIS  Work  Unit  No.  Civis  32221,  entitled  Probabilistic  Methods  in  Soil 
Mechanics.  Mr.  Richard  Davidson  was  the  OCE  Technical  Monitor. 

The  report  presents  an  introduction  to  statistical  quality  control  as 
applied  to  the  construction  inspection  of  engineered  embankments.  It  is 
intended  to  be  introduction  to  potential  users  who  have  little  or  no  background 
in  statistics.  Examples  in  the  report  are  drawn  from  actual  construction 
records  of  dam  projects,  and  IBM-compatible  microcomputer  software  supporting 
this  report  has  been  developed  under  separate  funding.  TVo  other  instructional 
reports  were  prepared  under  the  same  contract,  "Statistical  Analysis  of 
Geotechnical  Data,"  and  "Errm  Analysis  for  Geotechnical  Engineering,"  in 
addition  to  a  final  report. 

Ms.  Mary  Ellen  Hynes-Grif f in.  Earthquake  Engineering  and  Geophysics 
Division  (EEGD),  Geotechnical  Laboratory  (GL),  WES,  was  the  Contracting 
Officer's  Representative  and  WES  Principal  Investigator  for  CWIS  Work  Unit 
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STATISTICAL  QUALITY  CONTROL  OF  ENGINEERED  EMBANKMENTS 

PART  I :  INTRODUCTION 

Background 

Concern  with  quality  and  the  control  of  manufacturing  or  construction 
processes  to  assure  quality  underlie  modern  engineering  and  production. 

Indeed,  quality  control  is  as  old  as  engineering  itself.  On  the  other  hand, 
s  tatistical  quality  control  is  relatively  recent.  In  the  United  States, 
statistical  quality  control  first  came  into  its  own  with  the  wartime  production 
effort  of  1939-1945,  the  main  impetus  of  this  push  having  been  Army  Ordinance 
and  the  War  Production  Board.  The  military  influence  has  been  important  to  the 
introduction  of  statistical  quality  control  to  American  industry. 

Quality  control  in  construction  has  characteristics  which  are  both  similar 
to  and  different  from  quality  control  in  manufacturing.  Tbe  control  of  quality 
in  dam  projects,  especially  concerning  the  placement  and  compaction  of 
engineered  embankments,  is  critical  to  the  safety  and  performance  of  the  entire 
project.  Consequently,  a  well  planned  inspection  program  is  considered 
essential  on  any  moderately  large  project.  Current  CE  guidance  for  quality 
control  of  engineered  embankments  is  contained  in  EM  1110-2-1911,  dated  17 
January  1977.  It  is  not  s ta t i s t ica 1 ly  based,  but  is  experience  based. 

Purpose 

The  purpose  of  this  report  is  to  provide  potential  users  of  statistical 
quality  control  of  engineered  f l  1  Is  with  an  introduction  to  practical  concepts, 
definitions,  and  techniques.  'Itie  report  presents  simple  techniques  which  are 
intended  f  >r  use  by  readers  having  limited  familiarity  with  statistical  theory. 


The  report  does  not  attempt  to  survey  the  literature  of  statistical  quality 
control,  but  concentrates  on  a  few  chosen  techniques  that  fill  the  needs  of 
geotechnical  engineering  practice. 

General  Description  of  Statistical  Quality  Control 

The  placement  of  compacted  fills,  like  any  manufacturing  or  construction 
process,  varies  with  time.  The  physical  properties  of  soils  being  placed 
varies  in  moisture  content,  gradation,  plasticity  and  other  ways;  and  the 
process  of  placing  soils  varies,  for  example,  in  lift  thickness,  compactive 
effort,  and  climatic  conditions.  These  variations  cause  physical  properties  of 
a  resulting  fill  to  differ  from  one  point  to  another.  A  field  inspection 
program  intends  to  ensure  that — to  an  acceptible  level  of  conf idence--the 
completed  fill  conforms  to  specified  standards  and  thus  will  perform  its 
function  acceptably. 

Ideally,  an  inspection  program  could  non-destructively  screen  all  soil 
placed  in  a  fill  and  reject  those  materials  with  engineering  properties  not 
conforming  to  specified  standards  of  strength,  deformability  or  permeability. 
Such  a  program  would  guarantee  perfection  by  detecting  all  parts  of  a  fill 
which  were  flawed.  Unfortunately,  cost  and  the  lack  of  reliable  testing 
technology  precludes  this  ideal  program.  Instead,  typical  inspection  programs 
consist  of  limited  numbers  of  small-scale  tests  spread  thinly  throughout  a 
fill.  The  properties  measured  by  most  of  these  tests — for  example,  moisture 
content  and  dry  density — are  merely  surrogates  for  the  engineering  properties 
of  actual  interest,  although  some  engineering  properties  are  also  measured 


Statistical  quality  control  uses  simple  probalility  theory  to  develop 
inspection  sampling  plans.  These  plans  make  efficient  use  of  resources,  and 
can  be  related  to  a  quantitative  confidence  in  the  quality  of  a  finished 
product. 

Organization  of  This  Report 

This  report  is  organized  in  six  parts.  After  the  introduction,  Part  II 
summarizes  fundamentals  of  probability  and  statistics  which  are  necessary  for 
later  presentations.  Part  III  presents  basic  concepts  of  statistical  quality 
control  including  sampling  theory.  Part  IV  briefly  reviews  field  control  of 
compaction  operations.  Part  V  presents  quality  control  chart  techniques. 
Finally,  Part  VI  discusses  the  design  (i.e.,  planning)  of  sampling  schemes  for 
field  use.  Following  each  chapter  are  tables  and  figures,  and  plates 
presenting  example  calculations. 


PART  II:  FUNDAMENTALS 


This  section  briefly  reviews  mathematical  concepts  underlying  statistical 
quality  control. 


Probability  Theory 

Probability  theory  is  a  branch  of  pure  mathematics.  It  is  logical  and 
internally  consistent  in  the  sense  that  all  the  mathematics  of  probability 
theory  can  be  derived  from  a  small  set  of  axioms.  In  essence,  the  axioms 
specify  properties  that  "probability"  must  have,  for  example  probability  is  a 
real  number  between  zero  and  one.  Yet,  nowehere  do  the  axioms  say  what  the 
concept  of  probability  means.  As  a  result  many  interpretations  of  what 
probability  means  are  in  common  use. 

Frequency 

In  statistical  quality  control,  probability  is  usually  interpreted  to  be 
the  frequency  of  occurrence  of  some  event  in  a  long  series  of  similar  trials. 

A  trial  is  an  individual  occurence  producing  an  outcome  of  some  sort.  For 
example,  each  individual  lift  of  soil  placed  in  a  compacted  embankment  might  be 
considered  a  trial.  The  frequency  of  soils,  having  low  moisture  content  among 
these  lifts  (i.e.,  among  the  trials)  would  be  the  probability  of  soil  with  low 
moisture  content. 

Subjective  Probability 

An  alternative  interpretation,  common  in  geotechnical  engineering,  holds 
that  probability  is  a  rational  degree  of  belief.  Itie  probability  that  an 


important  solution  cavity  exists  in  a  limestone  dam  abutment  is  typical  of 
geotechnical  problems  which  cannot  be  easily  approached  using  the  frequency 
definition  of  probability.  Such  probabilities  have  to  do  with 
one-time  events,  past  experience,  and  amounts  of  information.  They  are 
personal  and  subjective.  They  are  not  related  to  frequencies,  actual  or 
conceptual. 

In  this  report  the  frequency  definition  of  probability  is  used,  for  it  is 
appropriate  to  quality  control  problems. 

Randomness 

A  key  concept  of  the  frequency  approach  to  probability  is  randomness. 

There  are  two  places  where  the  concept  of  randomness  is  important.  One  is  the 
description  of  a  construction  process  as  operating  in  a  random  manner;  the 
other  is  the  design  of  a  random  sampling  plan. 

A  process  is  operating  in  a  random  manner  when  any  part  of  the  output  may 
be  viewed  as  typical  of  the  output  as  a  whole.  That  is,  when  perturbations 
show  no  discernable  pattern.  Usually  it  is  not  possible  to  demonstrate  that  a 
process  is  operating  randomly;  rather,  it  is  only  possible  to  do  the  reverse, 
to  demonstrate  that  a  process  is  not  random.  This  is  done  by  showing  that  the 
output  of  the  process  does  in  fact  have  a  pattern  to  it. 

A  process  that  is  operating  in  a  random  manner  has  elements  or  events  with 
definite  probabilites  of  occurrence.  These  probabilities  may  not  be  known,  or 
may  be  known  only  to  the  extent  that  data  are  available  from  which  to  draw 
estimates.  Because  the  elements  or  events  have  associated  probabilities, 
statistical  theory  and  methods  can  be  used  to  characterize  them. 


The  other  place  randomness  is  important  is  in  the  design  of  random 


sampling  plans.  A  random  sampling  plan  is  one  in  which  sampled  elements  are 
chosen  with  definite  probabilities,  but  without  a  predictable  pattern.  In 
large  samples  the  relative  frequency  with  which  elements  are  sampled  should 
approach  those  probabilities;  however,  the  collection  of  elements  which  make  up 
any  specific  sample  reflect  a  chance  distribution. 

Conditional  Probability  and  Independence 

In  quality  control,  probability  is  commonly  defined  as  the  relative 
frequency  with  which  a  certain  event  occurs  in  a  long  series  of  similar  trials. 
For  example,  if  there  are  N  elements  in  a  large  set,  of  which  na  share  a  common 
property  A,  then  the  probability  of  an  element  within  the  set  having  property  A 
is. 


If  some  of  the  elements  also  share  a  common  property  B,  and  if  the  number  of 
these  elements  is  nb,  then  the  probability  of  property  B  within  the  set  is, 


P(B) 
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Consider  now  that  some  elements  in  the  large  set  possess  both  property  A 
and  property  B.  Let  the  number  of  such  elements  be  nab.  Graphically,  the 
number  of  elements  possessing  property  A,  property  B,  or  both  can  be  depicted 


as  in  Fig.  1.  If  we  consider  only  those  n^  elements  having  property  B,  the 
fraction  of  these  also  having  property  A  is  n^/n^*  Itiis  relative  frequency  is 
called  the  conditional  probability  of  an  element  having  property  A,  given  that 
it  is  known  to  have  property  B.  The  conditional  probability  is  denoted, 


P ( A  B)  = 


By  analogy,  the  conditional  probability  of  property  B  given  property  A  would 


P(B|A)  =  -  -4- 

a 

The  event  that  an  element  in  the  population  possesses  property  A  is  said 
to  be  independent  of  the  event  that  the  element  possesses  property  B  when  the 
probability  of  A  is  unchanged  by  knowing  that  an  element  possesses  property  B. 
Mathematically,  A  is  independent  of  B  if 


P(A  B)  =  P ( A ) 


Knowing  that  the  element  possesses  property  B  in  no  way  influences  the 
probability  that  it  also  possess  property  A.  If  the  event  that  the  element 
possesses  property  A  is  also  independent  of  the  event  that  it  possesses 
property  B,  then  properties  A  and  B  are  said  to  be  mutually  independent. 


Multiplication  Theorem 


Two  theorems  of  probability  theory  are  basic  and  often  encountered  in 
statistical  quality  control.  These  have  to  do  with  the  relationships  among  the 
probabilities  of  distinct  events.  The  first  is  the  multiplication  theorem. 

The  multiplication  theorem  states  that  the  probability  of  two  mutually 
independent  events  occuring  simultaneously  is  the  product  of  their  individual 
probabilities.  If  elements  possessing  properties  A  and  B  within  some  large  set 
are  mutually  exclusive  events,  then  the  probability  of  an  element  possessing 
both  property  A  and  property  B  is, 

P ( A  and  B)  =  P ( A)  P ( B)  .  -6- 

If  A  and  B  are  not  mutually  independent,  the  more  general  form  of  the  multipli¬ 
cation  theory  states  that  the  probability  of  them  occuring  simultaneously 
depends  on  the  conditional  probabilities. 


P ( A  and  B)  =  P(A)  P<B|A) 
=  P( B)  P( a| B) 
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Addition  Theorem 

The  addition  theorem  states  that  the  probability  of  either  one  or  the 
other  of  two  events  A  and  B  occuring  equals  the  sun  of  the  individual 
probabilities  of  their  occurrence,  minus  the  probability  that  they  occur 
toge  ther . 


P ( A  or  B) 


P  (  A )  +  P  (  B) 


P ( A  and  B) 


-B- 


I«.  it. 


H.  |I.  |t. 


In  Fig.  1,  if  area  is  taken  to  represent  probability,  the  addition  theorem  can 
be  considered  simply  a  statement  of  geometry.  The  area  contained  by  the 
combination  of  events  A  and  B  equals  the  sian  of  their  individual  areas,  less 
one  times  the  area  of  their  overlap  (i.e.,  P(A  and  B)  )  which  would  otherwise 


be  double  counted. 


Frequency  Distributions 


The  variability  of  data  on  production  output,  soil  properties,  or  other 
variables  is  com  niently  summarized  in  a  frequency  distribution,  the 
fundamental  tool  used  by  statisticians. 


Discrete  and  Continuous  Variables 

Fig.  2  shows  the  variability  of  standard  penetration  test  blew  counts 
measured  in  40  borings  in  a  silty  sand  deposit  at  a  dam  site.  Blow  counts  can 
only  assume  interger  values,  and  therefore  are  said  to  be  discrete  variables. 
Fig.  3  shows  variability  of  water  content  measured  in  R-tests  on  73  specimens 
of  a  compacted  clay.  These  strength  data  may  assune  any  real  number  value 
within  a  broad  range,  and  are  therefore  said  to  be  continuous  variables. 
Quality  control  in  geotechnical  engineering  must  deal  with  both  discrete  and 
continuous  variables,  and  many  methods  of  statistical  quality  control  apply  to 
each  in  a  similar  way. 


Histograms  and  Frequency  Distributions 

A  convenient  way  to  graphically  represent  scattered  data  is  in  a 
histogram.  A  histogram  graphs  the  number  of  measurements  falling  within 
specific  intervals  of  value  as  a  vertical  bar.  Fig.  2  shows  a  histogram  of  SPT 
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data.  For  obvious  reasons,  a  histogram  is  sometimes  called  a  bar  chart.  The 
height  of  the  bar  above  each  interval  shows  the  number  of  measured  values 
within  the  interval,  and  the  sum  of  the  heights  of  the  bars  equals  the  total 
number  of  measurements.  Fig.  3  shows  a  histogram  of  R-test  data. 

The  histogram  of  Fig.  3  divides  the  data  into  intervals  of  1%.  The  choice 
of  intervals  is  arbitrary,  but  the  intervals  should  be  of  uniform  width  and 
have  convenient  end  points.  If  too  many  intervals  are  chosen  the  general 
picture  of  relative  frequencies  will  not  be  obtained,  while  conversely,  if  too 
few  intervals  are  chosen  the  general  picture  will  be  blurred.  A  common 
rule-of-thumb  is  to  use  about  10  intervals.  More  detailed 
discussion  is  presented  in  the  report  "Data  analysis  for  geotechnical 
engineering"  (January  1986). 

A  frequency  distribution  is  constructed  from  a  histogram  by  dividing  each 
vertical  bar  by  the  total  number  of  measurements.  This  gives  the  relative 
frequency  of  observed  value  in  each  interval  as  a  decimal  fraction.  The  sum  of 
the  heights  of  the  bars  in  a  frequency  distribution  is  1.0.  Fig.  4  shows  the 
frequency  distribution  (right  had  side  scale)  corresponding  to  the  histogram  of 
Fig.  2.  Fig.  5  shows  the  frequency  distribution  corresponding  to  the  histogram 
of  Fig.  3. 

Cummulative  Distribution 

A  cummulative  distribution  of  discrete  or  continuous  data  is  constructed 
by  summing  relative  frequencies  starting  at  the  lower-value  end  of  the  data 
and  proceeding  toward  the  upper  value  end.  The  cummulative  distribution 


denoted  F(x)  gives  the  fraction  of  measurements  less  than  or  equal  to  a 


particular  value, 

F(x)  =  fraction  of  measurements  <.  x.  -9- 

Cummulative  frequencies  for  the  data  of  Figs.  2  and  3  are  shown  in  Figs.  6 
and  7.  The  cummulative  distribution  has  the  properties  that, 


For 

x  =  lower 

limit 

(or  -®) 

F(x)  =  0 

For 

x  =  upper 

limi  t 

(or  +a> ) 

-► 

F(x)  =  1 

For  discrete  data  the  cummulative  distribution  is  a  step  function  increasing  to 
the  right.  For  continuous  data  the  cummulative  distribution  is  typically  a 
smooth  S-shaped  curve. 

Importance  of  Frequency  Distributions 

Frequency  distributions  give  a  summary  view  of  the  variation  in  a  set  of 
data.  The  shape  of  the  distribution  suggests  whether  the  data  have  any  central 
tendency,  and  if  so,  where  along  the  x-axis  the  data  are  concentrated.  The 
width  of  the  distribution  indicates  the  dispersion  or  scale  of  variation  of  the 
data. 

Some  frequency  distributions  have  one  point  of  concentration  and  are 
thus  called  unimodal.  Others  have  more  than  one  and  are  called  multimodal. 
Usually,  soils  data  are  unimodal.  Multimodal  distributions  may  indicate  a 
mixture  of  data  from  different  soil  types  or  different  construction  procedures, 
that  is,  nonhomogeneous  data. 


The  frequency  distribution  also  shows  whether  the  variation  in  data  is 
symmetric  or  asymmetric,  that  is,  whether  high  and  low  variations  are  evenly 
balanced.  For  data  that  are  asymmetrically  distributed,  large  variation  from 
the  central  tendency  of  the  data  set  are  more  frequent  on  one  side  of  the 
center  than  on  the  other.  Ihis  is  illustrated  in  Fig.  8. 

Summary  Statistics 

Frequency  distributions  are  convenient  representations  of  data  for  visual 
inspection,  but  often  numerical  measures  of  distribution  characteristics  are 
useful  for  calculation  or  for  setting  standards.  Numerical  measures  are 
essential  for  developing  quality  control  criteria  and  quality  control  charts. 
The  most  important  numerical  measures  pertain  to  the  central  tendency  of  data 
and  to  dispersion. 

The  term  "statistic"  refers  to  any  mathematical  function  of  a  set  of 
measured  data.  For  example,  given  the  measurements  x-j,...,  xu,  any  function  y 
=  T(x-|,...,xu)  is  said  to  be  a  statistic  of  the  data.  The  arithmetical  average 
is  such  a  function,  the  largest  value  xmax  or  the  smallest  value  xm£n  is  such  a 
function,  and  so  on.  Any  of  these  ways  of  summarizing  the  data  would  be  called 
a  statistic.  Obviously,  there  are  an  infinite  number  of  statistics  which  could 
be  calculated  from  any  set  of  data,  but  the  most  useful  have  either  to  do  with 
the  central  tendency  of  the  data  along  the  x-axis  or  to  the  dispersion  of  the 
data. 

Central  Tendency 

The  most  common  measures  of  central  tendency  are  the  mean,  median,  and 
mode.  Hie  mean  is  the  arithmetic  average  of  a  set  of  data.  The  median  is  the 
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value  for  which  half  the  observations  are  smaller  and  half  larger.  The  mode  is 


the  most  frequent  value  (Table  1). 

The  mean  of  a  set  of  n  data  x  =  {xj,  . ..,  xn} ,  denoted  mx,  is  defined  as 
the  arithmetical  average. 


The  mean  is  the  center  of  gravity  of  the  frequency  distribution  along  the 
x-axis,  as  shown  in  Pig.  9. 

The  median  of  the  set  of  data  x  =  {x^ ,  ...,  xn} ,  denoted  xg.5,  is  the 
value  of  xn  which  half  the  data  are  less  than  and  half  more  than.  The 
cummulative  distribution  evaluated  at  the  median  is  0.5, 

F(xo.5)  =0.5  -1  2- 

i 

)  The  median  is  the  midpoint  of  the  data,  when  listed  in  increasing  or  decreasing 

order.  Common  practice  in  the  case  of  an  even  number  of  data  is  to  define  the 
median  as  half  way  between  the  two  middle  data,  that  is,  those  of  rank  (n/2) 
and  ( n/2  +  1 ) . 

The  mode  of  the  set  of  data  jt  =  { x-j  ,  ...»  xn}  ,  denoted  xQ,  is  the  most 
often  observed  value.  This  is  the  value  of  x  having  the  highest  ordinate  on 
the  frequency  distribution. 


Dispersion 


The  most  common  measures  of  dispersion  are  the  standard  deviation,  range, 
and  inner  quar tiles. 

The  Standard  deviation  of  a  set  of  data  =  {x-|  ,  .  .  . ,  xn}  ,  denoted  sx,  is 

defined  as  the  root  mean  square  variability  of  the  data. 
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in  which  mx  =  the  mean  of  the  data.  The  denominator  (n-1 )  rather  than  (n)  is 
used  to  correct  a  statistical  bias.  In  estimatinq  the  standard  deviation  from 
data,  the  mean  is  usually  also  unknown.  Ttius,  the  mean  must  be  estimated  from 
the  same  data  as  the  standard  deviation.  TViis  causes  the  averaqe  squared 
variability  about  mx  to  be  smaller  than  it  should  be.  On  averaqe,  it  is  smaller 
by  a  factor  (n-l)/n.  Correcting  for  this  error  gives  Bqn.  13. 

The  coefficient  of  variation  of  a  set  of  data  is  the  standard  deviation 
divided  by  the  mean, 

x  =  %/mx  —  1  4- 

The  coefficent  of  variation  is  used  to  express  relative  dispersion. 

The  variance  of  a  set  of  data,  denoted  Vx,  is  the  square  of  the  standard 


devia  tion. 
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In  many  statistical  calculations  the  variance  is  a  more  convenient  term 
than  the  standard  deviation,  and  is  thus  widely  encountered  in  statistical 
quality  control.  The  variance  is  the  moment  of  inertia  of  the  frequency 
distribution  about  mx. 

The  range  of  a  set  of  data,  denoted  r,  is  the  difference  between  the 
largest  and  smallest  values, 


rx  -  xmax  “  xmin 


The  range  has  poor  statistical  properties  in  that  it  is  sensitive  to  extreme 
values  in  a  data  set,  however,  it  is  easily  evaluated  and  therefore  often 
useful . 

The  inner  guar  tiles  of  a  set  of  data,  denoted  xg  .  25  an(^  *0.75'  are  the 
data  values  for  which  one-quarter  of  the  data  are  smaller  and  one-quarter 
larger,  respectively.  The  quartiles  may  be  found  from  the  cummulative 
distribution  as 


F(xo.25)  =  0.25 
F  ( *0 . 75  >  =  0.75, 


The  interquartile  range,  denoted  rn  c , 


is  less  influenced  by  extreme  values  than  is  the  range  itself,  but  it  is 
correspondingly  more  troublesome  to  compute.  Various  summary  statistics 
applied  to  the  R-test  data  of  Fig.  3  are  evaluated  in  Plate  1. 

Association  Among  Uncertain  Variables 

When  dealing  with  two  or  more  soil  properties  the  uncertainties  in 
estimates  may  be  associated  with  one  another.  That  is,  the  uncertainty  in  one 
property  estimate  may  not  be  independent  of  the  uncertainty  in  the  other 
estimate.  Consider  the  problem  of  estimating  'cohesion'  and  'friction' 
parameters  of  a  Mohr-Coulomb  strength  envelope.  If  the  slope  of  the  envelope 
to  the  Mohr  circles  is  mistakenly  estimated  too  steeply,  then  for  the  line  to 
fit  the  data  the  intercept  will  be  too  low.  The  reverse  is  true  if  the  slope 
is  estimated  too  flat.  Ttius,  uncertainties  about  the  slope  and  about  the 
intercept  are  not  independent,  they  are  associated  with  one  another. 

The  correlation  coefficient  for  paired  data  =  { (x-| , y-j  ) , . . . ,  (xn,yn)}  is 
denoted  pXy»  and  defined  as, 


Pxy 
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-19- 


In  effect,  the  correlation  coefficient  is  equivalent  to  a  normalized  product 
moment  of  inertia  in  solid  mechanics.  It  expresses  the  degree  to  which  two 
parameters  vary  together.  The  correlation  coefficient  is  non-dimensional 


because  deviations  of  x  and  y  are  measured  in  the  same  units  as  their 


respective  means.  Ttie  value  of  pxy  may  vary  from  +1  to  -1.  Px/y=+1  implies  a 

strict  linear  relation  with  a  positive  slope;  px>y=-1  implies  a  strict  linear 
relation  with  a  negative  slope;  px,y=o  implies  no  association  at  all  (i.e., 
independence) . 

Quick  Estimates 

Often  one  wants  quick,  approximate  estimates  of  means,  standard 
deviations,  or  correlation  coefficients  from  limited  numbers  of  data.  Some 
shortcut  techniques  are  available  for  th'is  purpose.  These  provide  economies  of 
time  and  effort  while  causing  sometimes  only  minor  losses  of  accuracy  or 
precision. 

Shortcuts  for  Estimating  the  Mean 

Rather  than  using  Equation  11,  a  quick  and  often  good  estimate  of  the  mean 
can  be  obtained  from  the  median.  The  median  is  the  middle  value  of  a  data  set. 
For  example,  if,  say,  five  data  are  listed  in  ascending  order  x^  ,  X2,  X3,  X4, 

X5 ,  the  median  is  X3.  For  an  even  number  of  data,  say  n=6,  the  difference 
between  the  two  middle  data  is  halved  to  give  the  median,  that  is  (x3+X4>/2. 

For  data  scatter  which  is  symmetric  about  its  central  value  and  for  small 
numbers  of  data,  the  sample  median  is  a  good  estimate  of  the  true  mean.  On  the 
other  hand,  if  the  data  scatter  is  asymme tr ic--for  example,  if  there  are  many 
small  values  and  a  few  large  values--the  sample  median  is  not  such  a  good 
estimator  of  the  mean. 

A  second  shortcut  for  estimating  the  mean  is  taking  one-half  the  sum  of 
the  largest  and  smallest  measured  values,  (l/2)(xmflx  +  xm1n).  This  estimator 


is  sensitive  to  the  extreme  values  in  a  set  of  measurements,  and  thus 


fluctuates  considerably.  It  is  not  a  good  shortcut  estimator  and  should  only 
be  used  with  caution. 

Shortcuts  for  Estimating  the  Standard  Deviation 

Rather  than  using  Equation  1 3,  a  quicker  estimate  of  the  standard 
deviation  from  small  numbers  of  tests  can  be  made  from  the  sample  range 
rx=(  xmax-xmin)  •  The  ran9e  is  the  span  of  data  from  largest  to  smallest.  Like 
the  standard  deviation,  the  range  is  a  measure  of  dispersion  in  a  set  of  data. 
However,  the  relationship  between  the  standard  deviation  and  the  sample  range, 
on  average,  depends  on  how  many  tests  are  made.  To  obtain  a  best  estimate  of 
sx  from  the  range  of  data  rx  a  multiplier  Nn  is  used  which  depends  on  sample 
size  (Table  2).  The  best  estimate  of  the  standard  deviation  is  sx  «  Nnrx  (see 
Plate  2). 

As  for  the  sample  median,  the  range  is  a  good  estimator  of  the  standard 
deviation  for  small  n  and  symmetric  data  scatter.  Even  for  modest  n  it  remains 
fairly  good.  However,  for  asymmetric  data  scatter  the  range,  which  is  strongly 
affected  by  outliers,  is  not  a  good  estimator  of  the  standard  deviation. 
Fortunately,  with  the  notable  exception  of  hydraulic  parameters  such  as 
permeability,  most  geotechnical  data  display  symmetric  scatter.  In  the  case  of 
hydraulic  data  a  logarithmic  transformation  (Lee,  et  al.,  1983)  usually  makes 
the  data  scatter  symmetric,  and  again  the  median  and  range  become  convenient 


estimators. 


Shortcuts  for  Estimating  the  Correlation  Coefficient 
Calculation  of  correlation  coefficients  by  Bqn.  19  can  be  time  consuming 
and  tedious.  A  simple  and  quick  approximation  is  obtained  graphically  from  the 
shape  of  the  scatter  plot  of  y  vs.  x.  The  method  works  well  whenever  the 
outline  of  the  scatter  plot  is  approximately  ellipical,  and  works  even  with 
small  numbers  of  observations.  Using  Chatillon's  (1984)  term  and  prodedure, 
this  is  called  the  balloon  method: 


STEP  1:  Plot  a  scatter  diagram  of  y  vs.  x. 

STEP  2:  Draw  an  ellipse  (balloon)  surrounding  all  or  most  of 
the  points  on  the  plot. 

STEP  3:  Measure  the  vertical  height  of  the  ellipse  at  its 

center,  h,  and  the  vertical  height  of  the  ellipse  at  its 
extremes,  H. 

STEP  4:  Approximate  the  correlation  coefficient  as  r  a  /l  -(h/H)  . 

An  example  of  the  method  is  shown  in  Fig.  10.  Ihe  balloon  method  gives  a 
correlation  coefficient  of  0.81,  whereas  the  correlation  coefficient 
calculated  by  Eqn.  19  is  0.83.  Empirically,  the  method  works  well  for  r>0.5. 

Shilling  (1984)  suggests  a  similar  balloon  method  for  approximately 
estimating  the  correlation  coefficient: 

STEP  1:  Plot  a  scatter  diagram  of  (y-my)/sx  vs.  (x-mx)/sy. 

STEP  2:  Draw  an  ellipse  surrounding  all  or  most  of  the  points 
on  the  plot. 

STEP  3:  Measure  the  length  of  the  principal  axis  of  the  ellipse 
having  positive  slope,  D,  and  the  length  of  the 
principal  axis  of  the  ellipse  having  negative  slope,  d. 

STEP  4:  Approximate  the  correlation  coefficient  as  r  (  d2  -d^  )  /  ( 1)2  +d-  )  . 

This  methods  works  about  as  well  as  Chatillon's.  Fbr  the  data  of  Fig.  10 
Shilling's  method  gives  r  0.80. 
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Probability  Distribution 


For  many  problems  in  statistical  quality  control  it  is  convenient  to 
approximate  the  empirical  frequency  distribution  for  some  category  of  data  by  a 
mathematical  function.  Surprisingly,  a  comparatively  small  set  of  mathematical 
functions  can  be  used  to  fit  a  broad  range  of  frequency  distributions 
encountered  in  the  field.  By  far  the  most  important  of  these  is  the  Normal  or 
bell-shaped  distribution.  Among  other  useful  distributions  are  the  log  Normal, 
Exponential,  and  4-parameter  Beta  distribution,  although  many  others  exist. 

The  Normal  distribution  is  discussed  here,  while  parallel  properties  of  the 
other  forms  are  given  in  Table  3. 

The  Normal  distribution  is  represented  by  the  equation 


,,  ,  1  2s' 

f(x)  =  — —  e  x  -20- 

/2n  s 

x 


in  which  mx  =  the  mean  of  x  and  sx  =  the  standard  deviation.  The  distribution 
is  unimodal  at  mx  and  symmetric  (Fig.  11).  Ihe  cummulative  distribution  of  x 
using  the  Normal  equation  is  found  from  the  area  under  the  frequency 
distribution  up  to  x, 


F( x )  =  fx  f ( x )  dx.  -21- 

—  CO 

The  normal  distribution  is  defined  for  -«><x<;+®,  but  the  area  under  the 


distribution  beyond  3  to  4  standard  deviations  from  the  mean  is  neglible. 

The  area  under  the  Normal  distribution  expressed  as  a  function  of  the 
standardized  variable 


and  calculated  by  Bqn.  21  are  given  in  Table  4.  Benjamin  and  Cornell  (1970) 
give  examples.  Numerically,  these  areas  can  be  approximated  by  the  series 
expansion  (Abramowitz  and  Segun,  1964), 

F(z)  =  fN(x)  ( b-,  t  +  b2t2  +  b3t3  +  b4t4  +  135 15 )  +  e  -23- 

in  which, 

=  0.319381530  t  =  (1+px)“l 

b2  =  -0.356563782  p  =  0.2316419 

b3  =  1.781477937  |e|  <  7.5  x  10"8  -24- 

b4  =  -1 .821 255978 
t>5  =  1.330  2  7  4  4  29 

and  ffj(x)  is  the  ordinate  of  the  Normal  distribution  function  evaluated  at  x. 
The  series  expansion  is  generally  more  convenient  than  Table  4  for  use  with 
computers. 

If  a  construction  process  is  operating  in  a  random  manner,  and  if  good 
estimates  of  the  mean  and  standard  deviation  are  available,  and  if  the  fre¬ 
quency  of  data  are  observed  to  be  well  modelled  by  a  representable  distribu¬ 
tion,  then  forecasts  can  be  confidently  made  about  the  future  performance  of 


that  process 


Ttiis  is  the  basis  for  statistical  quality  control 


Por  example, 


if  the  process  is  observed  to  produce  Normally  distributed  output,  then  a  chart 
such  as  Fig.  12  can  be  constructed  which  shows  the  process  mean  and  envelopes 
_K5sx  about  the  mean.  As  long  as  the  process  continues  to  operate  in  a  random 
manner,  and  the  mean,  standard  deviation,  and  frequency  distributions  remain 
unchanged,  then  a  confident  forecast  can  be  made  that  99.7%  of  the  output 
measurements  to  be  made  in  the  future  will  lie  within  the  +3sx  bound  (Figure 
12).  This  forecast  of  99.7%  comes  from  Table  4.  Such  forecasts  are  considered 
in  greater  detail  in  Parts  V  and  VI. 


Measure 


Symbol 


Formula 


Comments 


Central  Tendency 

Mean 

mx 

1/n  l  Xi 

center  of  gravity 

Median 

*0 .5 

F(xo.5)=0.5 

middle  value 

Mode 

*o 

x0=max  f(xi> 
i 

most  frequent  value 

Dispersion 

Standard  Deviation 

sx 

/  V  ( x . -m  ) 

/  n-1  L  i  x 

root  mean  square 
variation 

Variance 

Vx 

2 

sx 

moment  of  inertia 
about  mx 

Range 

rx 

xmax-xmin 

Interquartile 

Range 

r0 .5 

*0  .75_x0 .25 

i 


Ratio  of  average  range  to  standard  deviation  for  samples  from  a 
Normal  frequency  distribution. 


n 

Multiplier  Nn 

n 

Multipl: 

2 

0.886 

1  2 

0.815 

3 

0.591 

1  3 

0.300 

4 

0.486 

14 

0.294 

5 

0.430 

15 

0.288 

6 

0.395 

16 

0.283 

7 

0 . 370 

17 

0.279 

8 

0.351 

18 

0.275 

9 

0.337 

19 

0.271 

0 

0.325 

20 

0.268 

1 

0.315 

E 


Table  4 


—  Cumulative  frequencies  of  th<>  NjrmaL  distribution 
(from  Benjamin  and  Cornell,  11170). 


Cumulative  Probabilities  of  the  Normal  Probability  Distribution*  (areas  under  the 
norma)  curve  from  — to  z) 


~^T] 

.06 

.07 

.5199 

.5239 

.5279 

.5596 

.5636 

.5675 

.5987 

.6026 

.6064 

.6368 

.6406 

.6443 

.6736 

.6772 

.6808 

.7088 

.7123 

.7157 

.7422 

.7454 

.7486 

.7734 

.7764 

.7794 

.8023 

.8051 

.8078 

.8289 

.8315 

.8340 

.8531 

.8554 

.8577 

.8749 

.8770 

.8790 

.8944 

.8962 

.8980 

.9115 

.9131 

.9147 

.9265 

.9279 

.9392 

.9394 

.9406 

CO 

J- 

CO 

.9505 

.9515 

.9525 

.9599 

.9608 

.9616 

.9678 

.9744 

.9750 

.9756 

9798 

.9803 

.9808 

.9842 

.9846 

.9850 

.9878 

.9881 

.9884 

.9906 

9909 

.9911 

9929 

.9931 

.9932 

.9946 

.9948 

.9949 

.9960 

.9961 

.9962 

.9970 

.9971 

.9972 

.9978 

.9979 

.9979 

.9984 

.9985 

.9985 

.9989 

.9989 

.9989 

.9992 

9992 

9992 

.9994 

.9994 

9995 

.9996 

.9997 

9996 

.9997 

9996 

9997 

SUBJECT:  Example  calculations  of  Summary  statistics  for  (R-test)  data 


PROBLEM: 

Calculate  summary  statistics  from  experimental  R-test  data  on 
soil  strength. 

SOLUTION: 

1 .  Measured  data: 

Measured  values  of  R-test  data  as  shown  in  Figure  3. 

2 .  Measures  of  central  tendency: 

mean  —  mx  =  (1/n)  Zx^  =  32.6% 

median  —  xq . 5  =  32.1% 

mode  —  xq  =  28.5% 

3 .  Measures  of  dispersion: 


standard  deviation  —  1  .2  _  _ 

s  =  - -  Z(x.  -  m  )  =  6.5% 

v  n-1  1  x 


variance  — 

vx  = 

sx2  =  42.8% 

range  — 

r 

(xmax  ~  xmin)  =  29% 

f ractiles  -- 

*0.25 

=  27.4% 

*0  .50 

=  32.1% 

*0 .75 

=  36% 

interquartile  range 

—  *0.5 

=  (*0 .75  "  *0 .25>  =  8<6% 

PLATE  2 


SUBJECT:  Shortcut  estimates  of  summary  statistics. 

PROBLEM : 

Estimate  summary  statistics  using  shortcut  methods  and  compare  to  accurate 
calculations. 

DATA: 

Measured  Strength  (kPa):  38,  51,  43,  39,  48,  45,  42,  45,  49. 

ESTIMATE  MEAN: 

Shortcut  Method  Using  Median 

mx  <»  median  of  x, 

=  45  kPa 

By  Equation  2 

mx  =  -  I  Xi  =  ~  (400  kPa )  =  44.4  kPa 

x  n  x  g 

ESTIMATE  STANDARD  DEVIATION: 

Shortcut  Method  Using  Range 

w  =  ^xmax  ~  xmin) 

=  51-38  kPa 

=  1 3  kPa 

N  from  Table  1  (for  n=9) :  0.337 

sx  =  (0 .337)  (  1  3) 

=  4.4  k Pa 

By  Equation  3 

1 

sx  =  -  1  (Xi-uij,)'  =  4.2  kPa 


NUMBER  OF  DATA 


NUMBER  OF  DATA 


NUMBER  OF  DATA 


rnn 


e 

0 


Cumulative 

Frequency 


NUMBER  OF  DATA 


TOTAL  NUMBER  OF  DATA  •  225 

MEAN  •  S.9B2I  STD  DEv  •  4  4283 

number  of  outliers  beyond  lOuer  limit  of  plot 
number  of  outliers  beyond  upper  limit  of  plot 


N  VALUES  -  UNCORRECTED 

RE -REG  DAM  STA  25-00  -  STA  32-00 


Kiaure  '  —  Summary  s 


ry  statistics  of  a  frequency  distribution. 


,v»> 


Normal  or  bell-shaped  frequency  d 
(from  Chernoff  and  Moses,  1950,  E 
Decision  Theory,  John  Wl'ley  and  S 


MEASURED 


Quality  Assurance  and  Quality  Control 


The  terms  quality  assurance  ( QA )  and  quality  control  (QC)  are  used  in 
special  and  differing  ways  by  different  organizations.  In  this  report. 

Quality  assurance  means  an  inspection  program  aimed  at  assuring  that  soils 
placed  in  a  fill  meet  specifications. 

Quality  control  means  an  inspection  program  aimed  at  monitoring 
construction  performance  to  give  early  warning  of  changes  that  affect 
quality  and  thus  to  provide  a  basis  for  controlling  the  process. 

Quality  assurance  programs  prescribe  a  procedure  which  when  consistently 
applied  to  inspection  data  yield  a  specified  risk  of  accepting  lifts  of  given 
quality.  A  QA  program  provides  a  decision  procedure.  Quality  control,  on  the 
other  hand,  provides  a  way  of  estimating  lift  properties  and  the  changes  in 
those  properties  with  time.  A  QC  program  provides  a  monitoring  scheme.  QA 
provides  a  rule  by  which  the  owner's  risk  of  accepting  poor  quality 
construction  is  guaranteed  and  balanced  against  the  contractor's  risk  of  having 
good  quality  construction  rejected.  QC  provides  a  tool  by  which  owner  and 
contractor  alike  can  make  efforts  to  maintain  a  uniformly  high  quality 
product. 

Sampling 

Measurements  are  made  on  a  set  of  soil  specimens  or  at  a  set  of  locations 
in  order  to  estimate  the  properties  of  a  soil  deposit  or  an  engineered 
structure  as  a  whole.  Statisticians  call  this  set  of  measurements  a  sample. 


An  individual  piece  of  soil  is  called  a  specimen  to  distinguish  it  from  the 


concept  of  a  statistical  sample. 

The  soil  deposit  or  structure  whose  properties  are  of  interest  is  called 
the  target  population  (Fig.  13).  A  population  in  statistics  is  simply  a  larqe 
(or  infinite)  collection  of  elements.  Not  all  of  the  elements  in  the  tarqet 
population  may  be  accessible  for  samplinq.  Ttiose  that  are  accessible  are  said 
to  compose  the  sampled  population.  Prom  this  sampled  population  a  finite 
number  of  elements  are  selected  for  testing  and  this  set  is  called  a  sample. 

If  the  way  this  sample  is  chosen  satisfies  certain  rules,  the  sample  is  said  to 
be  a  probability  sample.  Statistical  methods  can  the  be  used  to  quantify  the 
uncertainty  in  estimates  from  the  sample  about  properties  of  the  sampled  \ 

population.  Statistics  is  powerless  to  say  anything  about  the  correspondence 

t 

between  sampled  and  target  populations,  however,  as  this  is  a  geological  or 
engineering  question.  ■' 

Scientific  Sampling 

The  concept  of  scientific  sampling,  or  probability  sampling,  is  central  to  [■ 

quality  assurance  and  control.  A  scientific  sample  is  planned  according  to 
statistical  principles.  The  importance  of  scientific  sampling  is  that  it  r 

allows  quantitative  statements  about  the  uncertainty  in  parameter  estimates 
which  result  from  sampling.  Other  samplinq  schemes--as  for  example,  ■ 

instructing  an  inspector  to  purposely  seek  out  areas  in  a  fill  that  appear 
poorly  compacted--certainly  have  merit  in  special  circumstances,  but  they  do  Q*, 

s 

not  alLow  the  quantitative  analysis  which  has  come  to  underlie  modern  v" 

V 

engineering  practice.  s‘ 

I 


To  be  a  probability  sample,  three  criteria  must  be  satisfied:  (1)  sample 


points  must  be  chosen  randomly,  (2)  all  elements  in  the  sampled  population  must 
have  a  non-zero  chance  of  being  sampled,  and  (3)  different  probabilities  of 
each  element  being  sampled  must  be  compensated  by  weighting. 

If  these  two  criteria  are  satisfied — and  only  if  they  are  satisfied — 
statistical  methods  can  be  used  to  determine  uncer tainties  properly  associated 
with  parameter  estimates.  This  means  that  for  statistical  methods  to  be  used, 
some  form  of  random  sampling  is  necessary.  Purposive  sampling,  by  which  an 
inspector  consciously  selects  for  testing  those  elements  that  appear  of  poor 
quality,  is  intuitively  appealing  and  can  provide  important  information,  but  it 
cannot  form  the  basis  for  statistical  quality  control.  From  a  purposive  sample 
there  is,  (a)  no  way  to  assign  quantitative  confidences  to  estimates  of  soil 
parameters,  (b)  no  way  to  explicitly  review  an  inspection  program  after  the 
fact,  and  (c)  no  way  to  establish  a  defendable  level  of  quality  assurance. 

People  also  talk  about  having  an  inspector  seek  out  a  'representative' 
sample.  Ibis,  too,  may  have  merit  in  special  circumstances,  but  it  does  not 
produce  a  sample  from  which  quantitative  conclusions  can  be  drawn.  No 
individual  sample  is  representative  of  a  sampled  population.  A  sample  contains 
specific  measurements  which  can  never  precisely  mirror  the  subtlety  of 
variations  in  the  sampled  population.  On  the  other  hand,  a  sampling  plan  can 
be  made  representative,  if  designed  by  scientific  principles,  in  that  it 
affords  every  element  within  the  sampled  population  an  equal  chance  to 
influence  estimates  that  are  made. 


Random  Sampling 


Scientific  sampling  requires  that  every  element  in  the  sampled  population 
have  a  non-zero  chance  of  appearing  in  the  sample.  It  does  not  require  these 
chances  to  aix  be  the  same,  only  that  the  relative  probabilities  are  known. 
This  condition  requires  that  elements  be  selected  from  the  sampled  population 
in  a  random  way.  Lacking  a  random  procedure,  the  assumption  that  each  element 
has  a  non-zero  chance  of  being  sampled  cannot  be  made,  and  the  relative 
probabilities  of  different  elements  being  sampled  cannot  be  assessed.  The  use 
of  a  random  procedure  attempts  to  avoid  any  form  of  association  between  the 
selection  of  elements  for  the  sample  and  the  properties  of  the  elements  that 
are  being  sampled.  Such  association  is  called  a  bias. 

Randomization  means  selecting  elements  of  a  sample  in  such  a  way  that  the 
two  conditions  of  probability  sampling  are  satisfied.  Randomization  can  be 
accomplished  many  ways.  A  conceptually  simple  but  operationally  clumsy  way  is 
to  pick  sampling  locations  by  a  random  number  generator  or  table  of  random 
numbers  (Table  5).  If  performed  faithfully  this  scheme  gives  each  element  in 
the  sampled  population  an  equal  chance  of  being  sampled.  Another  way  to 
provide  randomiza tion  is  to  layout  measurements  on  a  fixed  grid  and  then 
randomly  locate  the  first  point. 

Sampling  Plans 

An  essentially  infinite  number  of  sampling  plans  for  quality  assurance  or 
quality  control  satisfy  the  properties  of  probability  sampling.  These  are  all 
randomized  sampling  plans  in  the  sense  that  the  exact  elements  which  are 
sampled  depend  on  the  outcome  of  some  chance  event. 


A  more  convenient  sampling  plan  is  to  layout  sample  points  on  a  grid,  and 


then  locate  the  grid  in  the  field  by  randomly  selecting  its  first  point  (Fig. 
14b).  Only  one  pair  of  random  numbers  needs  to  be  chosen  from  which  all  of  the 
sample  points  are  determined.  The  disadvantage  of  a  grid  plan  compared  to 
purely  random  plans  is  that  any  spatial  periodicity  in  the  compaction  process 
may  bias  the  outcome.  An  advantage  compared  to  the  purely  rand'im  plan, 
especially  with  small  sample  sizes,  is  that  uniform  coverage  of  the  site  is 
assured . 

To  provide  coverage  while  at  the  ne  time  limiting  the  possible  effects 
of  periodicities,  stratified  random  sampling  plans  are  sometimes  tisel.  Using  a 
stratified  plan  the  sampled  area  is  first  divided  into  a  regular  array  of 
squares  or  rectangles  (Fig.  14c)  and  then  a  sample  p>int  is  randomly  located  m 
each. 


Another  common  plan  is  nested  sampling.  Nested  sampling  uses  a  pre-fixed 
grid  of  sample  points  with  varying  spacmgs  (Fig.  141).  The  first  punt  is 
located  randomly  as  in  grid  sampling  anl  fr'xn  that  punt  ill  the  res*  ire 
specified.  'Hie  principal  use  of  nested  sampling  is  for  estimating  sie  uil 
aspects  of  the  spatial  structure  it  soils  lata,  name  ly  the  nitocor re  la ♦  ion 
function  or  variograin,  (see  the  rep>r  t,  "Statistical  inalyu-  f  je  *e-iuu  « 
data".  Instructional  Repirt  wI.-ST).  The  u  ;e  >f  nested  sampling  in  pin’,  ity 
control  or  quality  assurance  of  compacted  fills  is  m  is  *  1  y  t  .r  toe  j  u  pis.-  >f 


assessing  me  asm  >*me  n  t  errors  >  r  nose  in  t  n»*  r  s.i  1  *  1  •:  a  i  i  *  i . 

Rand'im  •  •  1  unpe  l  sampl  l  r,u  plat,  s  i  r  d  t  -  n  u  se  «  ve  -  >  ,  1  ir  «*•  . .  i  *  i  t  ! 

extent  mu-;*  :>e  sanpl  e  1,  i  w  ,•  •  *  f  m  ;  •  i  i  i  .  r- ; 


operation  is  Large  compared  to  the  incremental  cost  of  testing.  Clumped 
sampling  involves  two  stages.  In  the  first  stage  a  number  of  seed  points  are 
randomly  chosen.  In  the  second  stage  a  number  of  sampling  points  are  chosen  in 
the  vicinity  of  each  seed  point.  At  both  the  first  and  second  stage  the 
sampling  plan  can  be  purely  random,  gridded,  stratified,  nested,  or  so  forth. 


The  sampling  plans  reviewed  here  are  typical  of  the  very  large  number  of 
possible  sampling  places.  In  practical  situations  the  constraints  of  a  parti¬ 
cular  project  may  dictate  that  a  specialized  plan  be  developed.  This  is 
accepted  practice  as  long  as  the  principles  of  probability  sampling  are  adhered 
to.  These  principles  dictate  three  things,  (1)  that  sample  points  be  chosen 
according  to  some  random  process  and  not  be  affected  by  the  intuition  of  an 
insjiector,  (2)  that  aLL  elements  within  the  population  to  be  sampled  have  a 
non-zero  chance  of  being  sampled,  and  (3)  that  if  the  probabilities  of  each 
element  being  s.impled  are  not  .ill  the  same,  these  differences  in  probability  be 
appropriately  compensated  for  by  weighting  when  the  da ta  are  analyzed. 

for  most  quality  control  and  quality  assurance  sampling  in  geotechnical 
engineering  the  probabilities  of  elements  within  the  s.impled  population  being 
s.impled  are  all  the  same.  Therefore,  for  these  sampling  plans  the  problem  of 
weighting  s. tin  pi  e  out  •••imes  is  seldom  of  concern.  For  those  cases  where 
w»*  i  fht  iii'i  i  i  . . .  iry,  Cochran  (  1  964  )  provides  techniques  an  1  practical 
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a  single  test  may  yield  more  than  one  numerical  result — for  example,  water 
content,  dry  density,  plasticity  index,  and  so  forth — but  for  now  nothing  is 
lost  by  considering  only  a  single  scalar  outcome. 

If  another  sample  of  n  specifications  is  now  taken  from  the  same  lift  of 
soil,  however  with  the  specimens  taken  at  slightly  different  places,  another 
set  of  n  numerical  data  will  result.  Each  of  these  will  differ  somewhat  from 
their  counterparts  in  the  first  sample,  because  the  soil  itself  varies  from  one 
spot  to  another  and  because  there  are  a  number  of  instrument  or  operator 
effects  which  influence  test  results.  'ttiis  variation  in  numerical  results  from 
one  sample  to  another  is  called  sample  variation.  Statistical  techniques  allow 
such  sample  variation  to  be  predicted  and  dealt  with  in  a  quantitative  way. 

The  sample  mean  (Eqn.  11),  sample  standard  deviation  (Eqn.  13),  and  other 
summary  measures  calculated  from  the  test  results  xi , . . . ,xn  are  simply 
mathematical  function  of  the  data.  If  the  data  vary  from  one  sample  to 
another,  so  will  the  summary  measures. 

Sampling  Variability  of  the  Mean 

The  sample  mean  mx  is  calculated  by  Eqn.  11.  If  many  tests  are  made 
(i.e.,  if  n  is  large),  variations  in  one  test  result  within  an  averaqe  will  be 
offset  by  variations  in  others,  and  as  a  result  mx  should  be  fairly  close  to 
the  actual  mean  of  the  sampled  population  mx.  In  this  report  the  actual 
sampled  populations  mean  is  denoted  by  a  prime,  mx' ,  as  compared  to  the  sample 
mean  which  is  denoted  without  a  prime,  mx.  On  the  other  hand,  if  few  tests  are 
made  (i.e.,  n  is  small),  variation  in  test  results  will  not  have  as  much 
opportunity  to  average  out,  and  as  a  result  the  sample  mean  may  deviate 
considerably  from  mx' .  This  sampling  variability  is  the  critical  factor  in 
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deciding  how  many  tests  must  be  made  in  a  quality  control  or  quality  assurance 


program . 

If  the  standard  deviation  of  the  sampled  population  is  known  and  if  the 
individual  measurements  are  independent  of  one  another,  then  the  means  of 
individual  samples  each  of  size  n  will  vary  with  a  standard  deviation  of 

smx=  sx  /  -25- 

For  exampLe,  Fig.  15  shows  a  histogram  of  sample  means,  each  corresponding  to 
different  sample  of  size  n  =  5  taken  randomly  from  the  SPT  blow  count  data  in 
Fig.  2.  The  standard  deviation  of  the  sampled  set  of  data  is  4.4  bpf,  while 
the  standard  deviation  of  the  variability  of  the  sample  means  is  2  bpf  «=  4.4 
bpf/  n.  If  plotted  as  a  frequency  distribution,  the  variability  of  the  sample 
means  will  be  approximately  Normally  distributed,  almost  without  regard  to  the 
shape  of  the  frequency  distribution  of  the  sampled  population.  The 
approximation  to  the  Normal  distribution  becomes  better  as  n  becomes  larger. 

In  the  more  common  case  the  true  standard  deviation  of  the  sampled 
population,  sx' ,  is  not  known,  and  thus  the  sample  standard  deviation,  sx,  is 
used  in  equation  25  to  approximate  the  variability  of  mx  about  mx' .  Using  sx 
rather  than  sx'  underestimates  the  variability  in  mx,  however,  because  the 
estimate  of  sx  differs  somewhat  from  the  true  standard  deviation  sx' .  To 
overcome  this  limitation  a  standardized  mean  is  used. 


in  which  mx  and  sx  are  the  sample  mean  and  sample  standard  deviation,  and  mx 
is  the  true  mean  of  the  sampled  population.  If  the  value  t  is  estimated 
separately  from  a  large  number  of  samples,  each  sample  containing  the  same 
number  of  observations,  the  frequency  distribution  of  t  over  these  many 
separate  samples  will  have  a  standard  deviation  of  1.0  and  a  shape  known  as 
Student1 s-t  distribution.  The  Student' s-t  distribution  looks  much  like  a 
Normal  distribution,  but  with  thicker  tails  and  a  higher  mode.  That  is,  the 
Student  t  has  somewhat  more  of  what  statisticians  call  kurtosis  than  a  Normal 
distribution  does.  Areas  under  the  Student  curve  are  given  in  Table  6,  and  may 
be  approximated  by  series  expansions,  as  given  by  (Abramowitz  and  Segun  (1964). 
The  shape  of  the  Student' s-t  distribution  and  thus  the  areas  beneath  it  depend 
on  the  number  of  measurements  within  a  sample,  n.  This  enters  Table  6  as  the 
degrees-of-f reedom  parameter  u  =  n-1  . 

Sampling  Variability  of  the  Standard  Deviation 

Just  as  the  sample  mean  varies  from  one  sample  to  another,  so  do  other 
summary  measures  such  as  the  sample  standard  deviation  or  sample  variance. 
Unfortunately,  the  statistical  results  for  the  variability  if  the  sample 
standard  deviation  and  variance  are  not  as  simple  as  those  for  the  sample  mean. 
For  samples  taken  from  Normally  distributed  data,  the  sample  standard  deviation 
varies  approximately  with  a  standard  deviation  of. 


in  which  n 


the  sample  size  (Snedecor  and  Cochran,  1980) 


The  sample 


variance,  sx2,  varies  approximately  with  a  standard  deviation  of, 


s 

s 


X 


2/2 
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x  /  n-1 


-28- 


Similar  results  are  available  for  non-Normally  distributed  data,  but  are  more 
complicated.  Most  basic  statistics  textbooks  discuss  these  results  (e.g., 
Snedecor  and  Cochran,  1980). 

Sampling  Variability  of  the  Range 

Because  the  sample  range  rx  =|xmax  -  xmin|  is  easier  to  calculate  than  the 
standard  deviation,  it  is  often  preferred  as  a  measure  of  variablity  in 
programs  of  quality  control.  For  a  sample  taken  from  Normally  distributed  data 
the  frequency  distribution  of  the  relative  range 

wx  =  rx/sx  -29- 

across  samples  of  size  n  is  tabulated  in  Table  7. 

In-Control  vs.  Out-of -Control  Processes 

The  above  discussions  are  based  on  the  concept  of  a  construction  process 
operating  in  a  random  manner.  When  a  process  is  operating  in  a  random  manner 
any  part  of  its  output  may  be  viewed  as  typical  of  the  output  as  a  whole  and 
perturbations  in  the  process  show  no  discernable  pattern.  A  construction 
process  operating  in  a  random  manner  and  producing  few  sample  outcomes  which 
deviate  subs  tan tia 1 ly  from  its  average  output  is  said  to  be  "in-control."  Pi  g . 
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16  shows  a  chart  of  compaction  data  for  a  process  operating  in-control.  The 


variations  in  these  data  appear  to  behave  randomly  without  trend  or  pattern. 

Few  construction  processes  operate  in-control  for  significant  lengths  of 
time,  and  even  when  a  process  is  in-control  minor  deviations  from  randomness 
have  to  be  overlooked.  The  statistical  theory  of  quality  control  and  quality 
assurance  is  based  on  the  idea  of  randomness  in  process  output,  and  thus  every 
effort  should  be  made  to  assure  that  non-random  factors  are  not  present. 

Whether  a  process  is  in-control  also  depends  on  the  level  of  detail  with 
which  the  output  is  scrutinized.  Obviously,  variations  in  the  output  of  a 
construction  process  are  not  truly  random.  At  some  level  of  detail  are  all 
explainable  by  physical  arguments.  The  notions  of  randomness  and  a  process 
being  in-control  have  to  do  with  engineering  decisions  and  the  cost 
ef fective-ness  of  further  reduction  in  output  variability.  Attempts  are  made 
to  identify  and  eliminate  all  major  sources  of  variability,  and  what  remains 
and  is  not  cost  effective  to  further  reduce  is  operationally  handled  as  if  it 
were  random  variation.  As  long  as  our  statistical  models  can  be  successfully 
used  to  portray  this  residual  variability  and  to  characterize  uncertainties 
which  arise  from  it,  then  the  process  is  for  engineering  purposes 
"in-control. " 


When  the  variability  in  the  output  of  a  construction  process  deviates  from 
randomness,  that  is,  when  significant  trends  or  patterns  begin  to  appear  in  the 
output,  statistical  models  no  longer  adequately  capture  the  important  features 
of  the  variations  and  the  construction  process  is  said  to  be  out-of-control. 


The  principal,  use  of  the  control  charts  of  Section  5  is  to  obtain  early  warning 
that  a  process  is  going  out-of-control,  and  to  identify  steps  that  might  he 
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Table  5  —  Table  of  uniform  [0,1]  random  numbers 
(from  Cochran,  1977). 


10 

27 

53 

96 

23 

71 

50 

54 

36 

23 

54 

31 

04 

82 

98 

04 

14 

12 

15 

09 

26 

78 

25 

47 

47 

28 

41 

50 

61 

88 

64 

85 

27 

20 

18 

83 

36 

36 

05 

56 

39 

71 

65 

09 

62 

94 

76 

62 

11 

89 

34 

21 

42 

57 

02 

59 

19 

18 

97 

48 

80 

30 

03 

30 

98 

05 

24 

67 

70 

07 

84 

97 

50 

87 

46 

61 

81 

77 

23 

23 

82 

82 

11 

54 

08 

S3 

28 

70 

58 

96 

44 

07 

39 

55 

43 

42 

34 

43 

39 

28 

61 

15 

18 

13 

54 

16 

86 

20 

26 

88 

90 

74 

80 

55 

09 

14 

53 

90 

51 

17 

52 

01 

63 

01 

59 

91 

76 

21 

64 

64 

44 

91 

13 

32 

97 

75 

31 

62 

66 

54 

84 

80 

32 

75 

77 

56 

08 

25 

70 

29 

00 

97 

79 

08 

06 

37 

30 

28 

59 

85 

53 

56 

68 

53 

40 

01 

74 

39 

59 

73 

30 

19 

99 

85 

48 

36 

46 

18 

34 

94 

75 

20 

80 

27 

77 

78 

91 

69 

16 

00 

08 

43 

18 

73 

68 

67 

69 

61 

34 

25 

88 

98 

99 

60 

50 

65 

95 

79 

42 

94 

93 

62 

40 

89 

96 

43 

56 

47 

71 

66 

46 

76 

29 

67 

02 

04 

37 

59 

87 

21 

05 

02 

03 

24 

17 

47 

97 

81 

56 

51 

92 

34 

86 

01 

82 

55 

51 

33 

12 

91 

63 

62 

06 

34 

41 

94 

21 

78 

55 

09 

72 

76 

45 

16 

94 

29 

95 

81 

83 

83 

79 

88 

01 

97 

30 

78 

47 

23 

53 

90 

34 

41 

92 

45 

71 

09 

23 

70 

70 

07 

12 

38 

92 

79 

43 

14 

85 

11 

47 

23 

87 

68 

62 

15 

43 

53 

14 

36 

59 

25 

54 

47 

33 

70 

15 

59 

24 

48 

40 

35 

50 

03 

42 

99 

36 

47 

60 

92 

10 

77 

88 

59 

53 

11 

52 

66 

25 

69 

07 

04 

48 

68 

64 

71 

06 

61 

65 

70 

22 

12 

56 

88 

87 

59 

41 

65 

28 

04 

67 

53 

95 

79 

88 

37 

31 

50 

41 

06 

94 

76 

81 

83 

17 

16 

33 

02 

57 

45 

86 

67 

73 

43 

07 

34 

48 

44 

26 

87 

93 

29 

77 

09 

61 

67 

84 

06 

69 

44 

77 

75 

31 

54 

14 

13 

17 

48 

62 

11 

90 

60 

68 

12 

93 

64 

28 

46 

24 

79 

16 

76 

14 

60 

25 

51 

01 

28 

50 

16 

43 

36 

28 

97 

85 

58 

99 

67 

22 

52 

76 

23 

24 

70 

36 

54 

54 

59 

28 

61 

71 

96 

63 

29 

62 

66 

so 

02 

63 

45 

52 

38 

67 

63 

47 

54 

75 

83 

24 

78 

43 

20 

92 

63 

13 

47 

48 

45 

65 

58 

26 

51 

76 

96 

59 

38 

72 

86 

57 

45 

71 

46 

44 

67 

76 

14 

55 

44 

88 

01 

62 

12 

39 

65 

36 

63 

70 

77 

45 

85 

50 

51 

74 

13 

39 

35 

22 

30 

53 

36 

02 

95 

49 

34 

88 

73 

61 

73 

71 

98 

16 

04 

29 

18 

94 

51 

23 

76 

51 

94 

84 

86 

79 

93 

96 

38 

63 

08 

58 

25 

58 

94 

72 

20 

56 

20 

11 

72 

65 

71 

08 

86 

79 

57 

95 

13 

91 

97 

48 

72 

66 

48 

09 

71 

17 

24 

89 

75 

17 

26 

99 

76 

89 

37 

20 

70 

01 

77 

31 

61 

95 

46 

26 

97 

05 

73 

51 

53 

33 

18 

72 

87 

37 

48 

60 

82 

29 

81 

30 

IS 

39 

14 

48 

38 

75 

93 

29 

06 

87 

37 

78 

48 

45 

56 

00 

84 

47 

68 

08 

02 

80 

72 

14 

23 

98 

61 

67 

49 

08 

96 

21 

44 

78 

37 

06 

08 

43 

37 

21 

34 

17 

68 

14 

29 

09 

34 

04 

58 

43 

28 

06 

36 

10 

43 

67 

29 

70 

44 

38 

88 

39 

54 

90 

69 

59 

19 

51 

41 

47 

10 

25 

62 

91 

94 

14 

63 

19 

80 

06 

54 

18 

66 

67 

72 

77 

63 

48 

59 

40 

24 

13 

27 

05 

90 

35 

89 

95 

44 

43 

80 

69 

98 

61 

81 

31 

96 

82 

42 

88 

07 

10 

05 

77 

94 

30 

05 

39 

83 

71 

46 

30 

49 

70 

52 

85 

01 

50 

25 

27 

99 

41 

28 

63 

61 

62 

42 

29 

68 

96 

83 

23 

56 

87 

83 

07 

55 

07 

49 

52 

83 

51 

14 

80 

62 

80 

03 

42 

86 

97 

37 

44 

22 

85 

39 

52 

85 

13 

97 

05 

31 

03 

61 

75 

89 

11 

47 

1 1 

09 

18 

94 

06 

19 

84 

08 

31 

55 

58 

79 

26 

88 

86 

30 

01 

61 

16 

96 

94 

46 

68 

05 

14 

82 

00 

57 

25 

60 

59 

24 

98 

65 

63 

21 

28 

10 

99 

00 

27 

89 

17 

95 

88 

29 

01 

84 

02 

78 

43 

07 

41 

08 

34 

66 

39 

68 

95 

10 

96 

32 

84 

60 

15 

31 

76 

58 

30 

83 

64 

47 

56 

91 

29 

34 

10 

80 

21 

38 

84 

00 

95 

01 

31 

76 

07 

28 

37 

07 

61 

20 

26 

36 

31 

62 

31 

56 

34 

19 

09 

98 

40 

07 

17 

81 

24 

33 

45 

77 

58 

01 

31 

60 

10 

39 

50 

78 

13 

69 

36 

90 

78 

50 

05 

62 

46 

72 

60 

18 

77 

47 

21 

61 

88 

32 

12 

73 

73 

99 

12 

02 

39 

56 

03 

46 

10 

62 

98 

19 

41 

19 

42 

74 

39 

91 

09 

24 

23 

00 

62 

44 

73 

67 

34 

77 

87 

29 

25 

58 

84 

05 

87 

31 

06 

95 

90 

56 

35 

03 

09 

17 

16 

29 

56 

63 

11 

16 

36 

27 

03 

68 

69 

86 

95 

44 

79 

57 

92 

36 

59 

22 

45 

44 

84 

1 1 

80 

45 

67 

93 

82 

53 

58 

47 

70 

93 

37 

68 

53 

37 

31 

77 

79 

13 

57 

44 

55 

66 

12 

62 

11 

27 

80 

30 

21 

60 

49 

99 

57 

94 

82 

97 

74 

06 

56 

17 

18 

83 

99 

47 

99 

41 

96 

53 

78 

72 

56 

12 

80 

73 

16 

91 

15 

79 

74 

58 

86 

50 

60 

00 

25 

12 

45 

57 

09 

09 

43 

12 

74 

49 

14 

38 

78 

94 

49 

81 

78 

86 

72 

04 

95 

84 

95 

48 

46 

45 

14 

93 

87 

81 

40 

24 

62 

20 

42 

31 

75 

70 

16 

08 

24 

85 

81 

56 

39 

38 

71 

26 

35 

03 

71 

59 

60 

10 

39 

66 

08 

99 

55 

64 

57 

10 

92 

35 

36 

12 

96 

88 

57 

17 

91 

Table  6 


Percentage  points  (i.e.,  doubie  tail  areas)  of  the  Student-t 
distribution.  For  areas  under  a  single  tail,  divide  by  two. 
From  Duncan,  1974. 
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4  318 

13 

■128  -259  -394  -538 

•694 

•870 

1079 

1  350 

1-771 

2-160 

2  650 

3  012 

4  221 

14 

•128  -258  -393  -537 

•692 

•868 

1076 

1.345 

1-761 

2  145 

2  624 

2  977 

4  140 

15 

■128  258  -393  -536 

691 

•866 

1  074 

1  341 

1-753 

2-131 

2  602 

2  947 

4  073 

16 

•128  258  -392  -535 

•690 

•865 

1071 

1-337 

1  746 

2  120 

2-583 

2  92  1 

4  015 

17 

•128  -257  392  534 

689 

■863 

1  069 

1333 

1-740 

2110 

2  567 

2  898 

3  965 

18 

■127  257  -392  534 

•688 

•862 

1067 

1-330 

1-734 

2-101 

2-552 

2  878 

3  922 

19 

•127  257  391  -533 

•688 

•861 

1066 

1-328 

1-729 

2  093 

2  539 

2  861 

3  883 

20 

•127  257  391  533 

687 

•860 

1  064 

1  325 

1-725 

2  086 

2  528 

2  845 

3  850 

21 

127  257  -391  532 

686 

•859 

1  063 

1-323 

1-721 

2  080 

2  518 

2  831 

3  819 

22 

■127  256  390  -532 

•686 

•858 

1  061 

1  321 

1  717 

2  074 

2  508 

2  819 

3  792 

23 

127  256  -390  532 

■685 

■858 

1  060 

1319 

1  -714 

2-069 

2  500 

2  807 

3  767 

24 

•127  256  390  531 

•685 

■857 

1  059 

1  318 

1-711 

2  064 

2  492 

2-797 

3  745 

25 

•127  -256  390  531 

•684 

■856 

1  058 

1-316 

1  -708 

2  060 

2-485 

2  787 

3-725 

28 

■127  256  -390  -531 

■684 

•856 

1  058 

1  315 

1  -706 

2  056 

2  479 

2  779 

3  707 

27 

•127  -256  -389  -531 

684 

•855 

1  057 

1314 

1  -703 

2  052 

2  473 

2  771 

3  690 

28 

•127  256  -389  530 

683 

•855 

1  056 

1  313 

1  -701 

2  048 

2  467 

2  763 

3  674 

29 

•127  -256  -389  -530 

•683 

•854 

1  055 

1311 

I  699 

2  045 

2  462 

2-756 

3  659 

30 

•127  -256  389  -530 

■683 

•854 

1  055 

1  -310 

1  697 

2  042 

2-457 

2  750 

3  646 

40 

•126  255  388  529 

■68 1 

■851 

1  050 

1  303 

1  684 

2  021 

2  423 

2  704 

3  -551 

60 

■126  254  -387  -527 

•679 

■848 

1  046 

1  296 

1  671 

2  000 

2  390 

2  660 

3  460 

120 

•126  -254  -386  -526 

•677 

■845 

1  -041 

1  -289 

1  658 

1  980 

2  -358 

2  617 

3  373 

«o 

•126  -253  385  -524 

■674 

•842 

1  036 

1  282 

1  645 

1  960 

2  326 

2  576 

3-291 

Table  7 


Percentage  points  of  the  distribution  of  relative  range 
w  =r  /s  for  small  samples  from  Normal  distributions 
(from  Duncan,  1974)  . 


Mean  w 

or 

di 

H 

Probability  That  tv  Is  Less  than  or  Equal  to  Tabular  Entry 

0.001 

m 

m 

0.025 

0.050 

0.950 

0.975 

0.990 

0.995 

0.999 

2 

1.128 

0.8525 

HI 

0.01 

0.02 

0.09 

2.77 

3.17 

3.97 

4.65 

3 

1.693 

0.8884 

119 

Esa 

0.19 

0.43 

3.31 

3.68 

4.42 

4 

2.059 

0.8798 

mm 

tsa 

0.43 

0.76 

3.63 

3.98 

4.40 

4.69 

5.31 

S 

2.326 

0.8641 

tu 

to 

0.66 

1.03 

3.86 

4.20 

4.60 

4.89 

5.48 

6 

2.534 

0.8480 

0.54 

0.75 

0.87 

1.06 

wm 

4.03 

4.36 

4.76 

5.03 

5.62 

7 

2.704 

0.833 

0.69 

0.92 

1.05 

1.25 

Bfttfl 

4.17 

4.49 

4.88 

5.15 

5.73 

8 

2.347 

0.820 

0.83 

1.08 

1.20 

1.41 

1.60 

4.29 

4.61 

4.99 

5.26 

5.82 

9 

2.B70 

0.808 

0.96 

1.21 

1.34 

1.55 

1.74 

4.39 

4.70 

5.08 

5.34 

5.90 

10 

3.078 

0.797 

1.08 

1.33 

1.47 

1.67 

1.86 

4.47 

4.79 

5.16 

5.42 

5.97 

11 

3.173 

0.7B7 

1.20 

1.45 

1.58 

1.78 

1.97 

4.55 

4.86 

5.23 

5.49 

12 

3.258 

0.778 

1.30 

1.55 

1.68 

1.88 

4.62 

4.92 

5.29 

5.54 

6.09 
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Populations  of  interest  in  samplinq. 


Figure  15  —  Histogram  of  the  nearis  >t 

from  a  large  data  set.  1 1;  all, 
samp  1 ed . 
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Thus,  a  sampling  inspection  program  has  two  goals:  to  control  homogeneity 
and  to  control  average  properties.  Which,  if  either,  of  these  goals  is  more 
important  depends  on  the  specific  situation. 

Target  Properties  vs.  Sampled  Properties 
The  properties  of  greatest  interest  to  the  engineering  performance  of  an 
embankment  or  fill  are  strength,  deformabi li ty ,  and  permeability.  However, 
these  target  properties  are  cumbersome  or  expensive  to  measure  directly,  so 
other  more  easily  or  quickly  measured  properties  are  used  in  their  place.  By 
far  the  most  commonly  sampled  properties  in  construction  inspection  of  fills 
are  compaction  water  content  and  dry  density.  The  fact  that  these  are 
correlated  to  strength,  deformabili ty ,  and  permeability  makes  them  useful 
sur roga tes. 


Tests  for  Water  Content  and  Dry  Density 
A  variety  of  tests  are  available  for  measuring  water  content  and  dry 
density.  Water  content  can  be  measured  directly  by  oven  drying  a  specimen  and 
determining  change  in  weight.  Dry  density  can  be  measured  directly  by 
ascertaining  the  weight  and  volume  of  a  specimen,  as  for  example,  with  a  sand 
cone  density  test. 

Water  content  and  dry  density  can  also  be  measured  indirectly  using 
various  devices,  for  example  by  nuclear  gage.  These  indirect  tests  are 
typically  less  expensive  than  direct  testing  but  also  less  accurate.  In 
certain  cases  economies  can  be  gained  by  combining  a  small  number  of  direct 
tests  with  a  larger  number  of  indirect  tests. 

')"'irr  i  pt  i  ins  of  field  density  and  water  content  tests  are  given  in  Lam  be 
(It  i),  She  r  a  r  i ,  et  a  1 .  ( 1  9b  1)  ,  Engineer  Manual  EMI  1 10-2-1 91  1 ,  USBR  (1960),  and 
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AASHO  and  ASTN  standard  test  procedures. 

Compaction  Specifications 

Specifications  for  compaction  quality  are  most  often  placed  on  water 
content,  dry  density,  or  both.  Fhr  example,  water  content  of  the  fill  at  time 
of  compaction  might  be  specified  to  be  within  +2%  of  standard  or  modified 
Proctor  optimum.  Dry  density  might  be  specified  to  be  at  least  95%  of  standard 
or  modified  Proctor  optimum.  These  are  performance  specifications. 
Specifications  are  also  placed  on  compaction  equipment  and  procedure.  Fbr 
example,  a  specified  number  of  passes  with  equipment  of  specified  minimum 
capacity  may  be  required  in  addition  to  some  specified  water  content  range.  As 
an  example,  on  the  USAE  Carters  Dam  Project,  Georgia,  compaction  specifications 
for  the  imperious  core  required  as-placed  water  contents  to  be  +2%  standard  or 
modified  Proctor  optimum  and  the  fill  to  be  compacted  by  a  minimum  number  of 
passes  using  specified  equipment.  If  placed  materials  were  found  to  have  water 
contents  more  than  ±2%  from  optimum,  the  cost  fell  to  the  contractor  to 
moisten,  dryout,  or  remove  the  material.  If  placed  materials  were  within  water 
content  specifications  and  had  been  properly  compacted,  but  were  less  than  95% 
standard  or  modified  Proctor  optimum  density,  then  the  cost  fell  to  the  owner 
to  undertake  additional  compaction  or  to  remove  the  material.  These  are 
compliance  specifications. 
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PART  V:  QUALITY  CONTROL  CHARTS 

Quality  control  techniques  are  used  by  both  the  contractor  and  the  owner 
to  monitor  the  progress  of  construction,  and  thus  to  quickly  identify  changes 
in  soils  or  operational  procedures  before  these  changes  adversely  affect 
construction  quality. 

Quality  control  differs  from  acceptance  sampling  in  that  QC  has  the 
principal  purpose  of  identifying  changes  in  construction  materials  or 
procedures  before  those  changes  adversely  affect  construction  quality.  When  a 
change  is  detected,  efforts  are  made  to  find  assignable  causes  and  fix  them. 
Acceptance  sampling,  in  constrast,  has  the  principal  purpose  of  assuring  that 
soils  placed  in  a  fill  meet  sj>ee  i  f  ioa  1 1  oris.  Based  on  a-.-eptan-e  sampling 
results  soils  are  either  accepted  or  rejected  as  part  >f  a  mill ♦y  assurin-e 
program . 
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many  chance  causes  w!i  i  oh  produce  variation  in  soil  Vnsity,  ai'H  stir-  ronfenf , 


or  ot.ner  properties  ire  similar  to  the  many  forces  which  cause  a  t  isse  1  coin  t 
land  up  tun ids  or  tails.  Such  variations  follow  predi  -tatile  laws  of 
pr  oba!)  i  1  1  tv. 

>n  *  ;k*  other  haul,  other  variations  in  quality  are  lue  to  assi  triable 
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regularity.  Therefore,  when  the  process  is  out-of-control  an  effort  is  made 
to  find  assignable  causes. 

Suppose  that  samples  of  fixed  size  n  (=  number  of  tests)  are  taken  from 
lifts  being  placed  in  a  compacted  fill.  From  each  of  the  n  tests  a  measured 
value  of  some  soil  property  results.  From  these  n  values  certain  statistics 
are  calculated,  for  example  the  sample  mean  mx={ 1/n ) ,  standard  deviation 
sx=|/(1/n-i>  (n^-mx)  ,  or  range  rx=(  xmax_xmin^  •  Being  sample  results,  these 
statistics  will  be  subject  to  fluctuation  from  one  sample  of  n  to  another. 
However,  if  the  variations  are  due  to  chance  causes,  the  frequency 
distributions  of  the  sample  mean,  standard  deviation,  range,  or  other 
statistics  are  known  to  follow  the  regular  distributional  forms  discussed  in 
Part  It.  For  example,  the  sample  mean  mx  is  known  to  have  a  frequency 
1 1 str lbut ion  in  repeated  sampling  wh i oh  is  approximately  normal  (exactly  normal 
if  the  soil  properties  being  tested  ,ir»  themselves  Normally  distributed).  The 
aver ige  value  >t  the  sample  mean  mx  equals  the  real  mean  mx' ,  while  the 

Utah  lari  !<*vi  t »  1  h,  s,nx,  of  mx  equals  sx'/  n.  From  Table  4,  only  0.2%  of  the 

•  ample  meati  »  shoal  i  lie  outside  i  fsmx  interval  iliout  mx'  . 

No.,  i  tu  ,  i ,  'ha  *  , '  i  *  i  u  f  i  -s  it  i'ii'  i.tmpl  •>  la  t  a ,  is  fir  example,  the  sample 
. .  '  >  i"  |  ;  *  i  u  1  i '  1  .•  i  i  ’  i  m ,  ire  len  i  r  .■  1  he  r  .•  w ;  •  tv  n  1 1  i  prime.  '-’or 

•  i  * '  i  ",  •  .  '  ir  (o-  ■!.'  n.  -■»  i  •  i  •  *  t  -s  ■  t  ♦  h**  wh  1  ••  lamp  1  "  1  |»  'pal  » »  i  m,  as  f  >r 
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Control  Charts  and  Control  Limits 


A  control  chart  is  a  device  by  which  the  state  of  statistical  control 
(i.e.,  that  a  process  is  in-control)  is  operationally  defined.  It  is  used  to 
attain  control  in  a  new  process,  and  check  that  control  is  maintained  in  an 
ongoing  process. 

A  control  chart  is  constructed  by  plotting  values  of  mx,  sx,  rx  or  other 
sample  statistics  as  a  function  of  time  or  of  some  other  dimension  for  ordering 
sample  results  (e.g.,  lift  sequence  number).  The  sample  statistics  are  plotted 
against  the  vertical  axis,  time  or  other  dimension  against  the  horizontal  axis. 
A  horizontal  line  is  drawn  through  the  actual  mean  mx' ,  which  could  be  fixed  by 
specification  or  calculated  from  date.  Two  other  horizontal  lines  are  drawn, 
one  above  mx'  and  one  below  mx',  showing  limits  which  are  highly  likely  to 
contain  the  samp-le  results.  These  are  the  control  limits:  the  upper  control 
limit  (IJCL)  and  the  lower  control  limit  (LCL).  Fig.  16  shows  a  typical  control 
chart  for  individual  compaction  data. 

If  sample  values  are  plotted  foi  a  substantial  range  of  production  and 
time,  and  if  all  these  values  fall  within  the  interval  formed  by  the  UCL  and 
LCL,  and  if  tne  data  show  no  cycles  or  runs,  then  it  is  concluded  that  the 
construction  process  is  in-control  for  that  particular  attribute.  If  the  data 
do  not  conform  '  this  pattern,  then  the  conclusion  is  drawn  that  variability 


in  the  constructed  pioduct  is  not  .  xplainahle  by  chance  factors  al  >ne  and  an 
aasi  in. tnle  ailed  a)  is  sought. 

1  r  i.ii  a  .fati'd  1  -Mil1:!  view,  control  limits  are  related  to  *  he  testing  .f 
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random  variation  is  reasonable  and  accepted.  If  the  results  would  be 

improbable  based  on  the  hypothesis  of  random  variation,  then  that  hypothesis  is 

rejected.  The  choice  of  control  limits,  and  the  associated  probability  of 

their  being  exceeded,  is  arbitrary.  For  example,  the  probabililty  could  be  set 

at  0.01,  implying  control  limits  of  ±  1.65  standard  deviations;  or  at  0.002, 

implying  control  limits  of  ±  3.0  standard  deviations.  Narrowing  the  control 

limits  means  increasing  the  risk  that  the  hypothesis  of  random  variation  will 

be  rejected  when  actually  it  is  true.  For  example,  with  limits  set  at  ±  3 

deviations  there  is  a  chance  of  0.002  that  a  process  actually  in  control  will 

fixed  at  (mx'+-3sm).  while  a  lower  control  limit  (LCL)  is  fixed  at  (mx'-3sm) 

Usually  the  standard  deviation  of  mx,  that  is,  s  ,  is  estimated  from  the  data 

mx 

as  s  =sx//n.  Sometimes  sx  itself  is  specified  as  a  target  homogeneity. 
mx 

Fig  17  shows  an  m-chart  for  compaction  control  data  on  a  dam  project. 

To  control  current  production,  a  sample  of  size  n  is  taken  periodically 
from  material  placed  in  the  fill  and  the  average  mx  of  the  n  test  results  is 
plotted  on  the  control  chart.  If  all  the  mx  lie  within  the  UCL  and  LCL,  the 
construction  process  is  concluded  to  be  in-control.  If  any  mx  falls  outside 
either  the  UCL  or  LCL  the  process  is  deemed  to  be  out-of-control.  When  the 
deviation  outside  one  of  the  control  limits  is  adverse,  for  example,  when  mean 
compacted  density  falls  below  the  LCL,  specific  cause  for  the  variations  are 
looked  for  with  the  intent  of  improving  the  construction  process  and  thus  the 
product  of  that  process.  When  the  deviation  beyond  a  control  limit  occurs  on 
the  favorable  side,  for  example,  when  the  mean  compacted  density  exceeds  the 


UCL,  either  no  action  is  taken  or  the  causes  of  this  unusually  high  quality  are 
searched  for  in  order  to  learn  how  to  permanently  improve  quality. 

Probabilities  of  individual  sample  means  exceeding  the  UCL  or  LCL  can  be 
found  by  reference  to  Table  4.  For  soil  property  data  which  are  themselves 

I 

m  -  m 
x  x 

Normally  distributed  the  probabilities  from  Table  4  are  exact  for  - - — 

s 

x 

For  soil  property  data  that  are  not  Normally  distributed--presuming  that  the 
distributions  are  not  bizarre — the  probabilities  of  Table  4  are  still 
approximately  correct  even  for  sample  sizes  as  small  as  3. 


Control  Chart  for  Sample  Range  vy 


A  control  chart  on  the  sample  range, 


rx  -  xmax  ~  xmin  -30- 

shows  variation  in  the  range  as  a  function  of  time.  Hie  central  line  on  an  R- 
chart  is  fixed  at  the  empirical  average  range  in  past  production,  or  in  special 
circumstances  is  set  by  specification  on  acceptable  variability  of  the 
compacted  fill.  Hie  control  limits  are  usually  set  at  1  3sr,  in  which  sr  is 
the  standard  deviation  of  the  sample  range.  Both  the  average  range  and  the 
standard  deviation  of  the  range  can  also  be  related  to  the  standard  deviation 
of  the  soil  properties  being  samples  sx' .  Fig.  IB  shows  an  R-chart  for 


compaction  control. 


If  data  fall  inside  the  UCL  and  LCL  on  an  R-chart,  the  construction 


process  is  deemed  to  be  in-control  with  respect  to  homogeneity.  When  a  single 


data  point  falls  outside  the  UCL  or  LCL  the  process  is  deemed  to  be  out-of¬ 


control  with  respect  to  homogeneity.  In  the  latter  case  actions  are  taken  to 


find  assignable  causes.  A  sample  result  above  the  UCL  is  usually  considered 


adverse  and  efforts  should  be  made  to  find  out  the  cause  of  the  variability  and 


fix  it.  A  sample  result  below  the  LCL 1 s  usually  considered  favorable  and 


efforts  can  be  made  to  find  out  what  is  being  done  so  well  so  that  the 


construction  process  can  be  improved. 


Because  an  m-chart  and  an  R-chart  control  for  different  aspects  of 


quality,  a  process  may  be  in-control  on  one  but  out-of-control  on  the  other. 


An  m-chart  controls  for  the  mean  or  average  quality  of  the  compacted  fill.  An 


R-chart  controls  for  the  uniformity  with  which  compacted  materials 


are  being  placed.  Compacted  soils  may  be  on  average  sufficiently  dense,  but 


unacceptabely  heterogeneous.  On  the  other  hand,  the  soils  may  be  on  average 


sufficiently  uniform  but  unacceptably  loose. 


Cumulative  Reject  and  Related  Charts 


Unlike  most  industrial  applications  of  quality  control  charts,  construc¬ 


tion  involves  a  single  project  with  a  clearly  identified  beginning  and  end.  As 


a  result,  certain  quality  control  charts  are  very  useful  in  construction  even 


though  they  are  not  widely  used  in  the  factory.  One  of  these  is  the  cumulative 


reject  chart. 


✓ 


Cumulative  Reject  Chart 


The  cumulative  reject  chart  plots  the  cumulative  number-;  of  tests  htvmq 
results  outside  specified  limits,  against  time  test  sequence  number,  or  a 
similar  indicator  of  test  order.  Fig.  19  shows  cumulative  reject  iatu  f  ir  a 
compaction  inspection  program  on  the  imperious  core  of  a  rook  fill  dam.  The 
upper  figure  (a)  shows  cumulative  rejects  due  to  inadequate  densities.  Hie 
middle  (b)  and  bottom  (c)  figures  show  numerical  values  jf  water  content  arid 
dry  density,  respectively,  for  the  rejected  tests.  These  are  plotted  alone 
with  the  cumulative  reject  test  so  that  the  cause  of  rejection  or  my  ti  ers  l  in 
the  cause  can  be  readiLy  seen. 

In  Fig.  19  cumulative  reject  is  plotted  against  test  sequence  number.  \s 
a  result,  the  slope  of  the  curve  gives  the  rate  of  rejects  at  my  print  iurini 
the  project.  In  all,  1  1 75  insjiection  tests  wore  made  on  t  he  imperious  -  re, 
which  39  were  rejected  either  for  being  outside  • 2*  Fr  icter  ,pt  jmur  wi*-e 
content  or  for  having  iry  d-'nsity  less  than  9F  t  Procter  opt.  imun .  The  iufe  -f 
rejects  for  the  entire  protect  is  3H/117S  -  3%. 

While  the  average  rite  .>{  reject,  tests  was  b*  f  >r  the  entire  i  n  ; 

cote,  dm  l  :tq  early  phases  of  construction  the  rate  wm  much  hi  :*;.»?  mi  tit  ;  •  : 
later  pnases  t  ht-*  rate  was  mu-di  1  wit.  *\t  the  s*  ir’  >f  the  •*. 

reject  ion  rate  re  ache  i  i  n  l  ;h  >\  about  X\%  ,  j  r  id  n  i  .  y  *  i :  t  ;  • :  :  :  •  •  ■  i  •  ■ . 

2*.  neai  n  nt|  let  1  m.  'Pi  j  ■;  ara  iua  I  im»  Ca-  a  ly  ie-Tease  in  • 


control.  Itie  mean  quality  was  low  and  there  was  considerable  variability  in 
compaction  properties.  As  construction  proqressed  the  process  was  brought  i nto 
tighter  control.  The  mean  quality  was  better  maintained  and  the  variability 
reduced . 

The  cumulative  reject  chart  can  be  used  to  monitor  a  number  of  subtle 
changes  in  the  construction  process.  Fig.  20  shows  schema cical ly  the  effect  of 
a  major  'hange  in  the  construction  process,  for  example,  a  change  of 
contractor,  change  of  equipment,  or  change  of  borrow  material.  The  change 
causes  a  break  in  the  smooth  progression  of  the  learning  curve,  usually 
starting  another  learning  cycle. 

careful  insi^'-ti.in  of  Fig.  11  shows  two  such  breaks.  The  first  occurs  at 
iDour  te s t  number  oh.  'Pie  reject  rate  for  tests  1  to  50  is  about  8%.  From 
t’o'  n  '  *  ■  abou*  H  t  he  reject  rate  increases  sharply  to  about  30%.  In  fact, 
lue  •  *  he  earning  effect  file  rate  of  rejects  should  be  expected  to  decrease 
n  ‘  ; ",  -re.i  i  :  •  moruct;  oi  pr  js.  ^  retrosjwctive  analysis  shows  t  ha  t  a  * 

ii  .  •  ’  ■  ■  •  iw  ~a  '  ••  r  i  i  w  i  ?  ’hang**  } .  Jn»*  borr  iw  sour  •*>  was  is.-  :  r. 

r  f  '  i  ‘  i  • :  l  i  s  -  i  r  v  o  ;  ..  ;e  i  Iron  ‘here  -in.  he.  'a,i  se  *  he  •  nan  je 


Moving  Average,  Standard  Deviation  and  Range  Charts 


Convenient  accessories  to  the  cumulative  reject  chart  are  the  movie.  : 
average,  moving  standard  deviation  and  moving  range  charts.  These  provi  le 
smoothed  information  on  changes  in  construction  output  from  which  trer.ls  ran  n 
more  easily  identified. 

In  a  standard  m-chart  the  averages  of  samples  of  n  tests  are  plotted  as  a 
function  of  time  or  some  other  ordering  index.  A  different  set  of  n  tests  is 
used  for  each  point,  and  the  assumption  is  made  that  each  test  is  independent 
of  every  other.  Thus,  each  plotted  sample  mean  mx  is  also  independent  of  the 
sample  means  of  adjacent  to  it,  presuming  that  the  construction  process  is 
in-control  and  operating  in  a  random  manner  (Note:  in  practice  the  problem  of 
serial  correlation  in  the  construction  process  itself  sometimes  arises,  but 
such  autocorrelation  is  beyond  the  scope  of  the  present  report) .  The  use  of 
-.-.-harts  typically  presumes  that  many  data  are  being  collected  and  that 
;  r  j-m  :  t  ion  output  is  fairly  high. 

- j r  many  cases  in  construction  the  rate  of  testing  is  more  modest  or 
■  - ;  •  •  .  proceeds  more  slowly.  Often  a  considerable  time  is  required  to 
'  .  -i  : v  i  iuu  L  lifts  that  are  to  be  tested.  In  these  cases  a  moving 

' .  i  r  •  -  ay  :>•-■  more  convenient  than  the  standard  m-chart.  The  moving 

-  '  .  i  ■  •»  i  true  ted  in  the  same  way  as  an  m-chart,  but  it  provides  a 

"  « -  now  the  construction  process  may  be  changinq. 

ir-  cal ou i a  ted  over  windows  of  fixed  size  n.  For  a  long 
<•,...,  •<y )  ,  an  averaeg  is  calculated  for  the  first  n 
'  •  i •  r  a  is  calculated  for  the  next  n  data 


distribution,  as  illustrated  in  the  above  example,  may  indicate  a  change  in  the 
construction  process  which  needs  to  be  monitored  or  controlled. 

Cumulative  Sum  (CUSIJM)  Chart 

Changes  in  the  output  of  a  construction  process  are  sometimes  more  quickly 
detected  by  monitoring  the  change  from  one  test  to  the  next  rather  than  the 
absolute  value  of  the  test.  For  example  the  changed  conditions  which  appear  in 
Figs,  ig  and  21  become  apparent  earlier  when  increments  of  test  results  are 
p 1 ,:> t ted . 

The  most  oinmon  way  to  monitor  increments  is  by  the  cumulative  sum  or 
ousutn  chart.  Hie  cusum  -'hart  uses  trends  m  the  gc  data  to  identify  process 
changes,  rather  than  treating  the  data  from  each  lift  by  themselves.  Hie  maior 
advantage  of  ;  us  urn  charts  over  in-  or  r-eharts  (sometimes  Snewhar  t  charts  after 
tie*  sta  t  j  st.  ician  who  first,  proposed  them)  are  t.ha  f  'hanges  are  identified  more 
pil  :k  1  y  ,  par  t l  cu  1  a  r  ly  mo.lest  'hanges,  and  that  the  tine  or  location  of  the 

hanje  mu  be  precisely  determined. 

'V  "rein  •hir’  -,f  i  s.-gacn-'y  if  la  'a  x ,  .  .  .  ,  x  ,  ,  .  .  .  i  s  -ons f  ru«' fed  by 
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As  lor.  }  as  the  construction  process  is  m  control ,  sue 
should  lie  on  a  .straight  line.  Sone  individual  x,  will  lie 
in  which  case  the  cusuni  curve  will  rise;  and  som*1  i  n  i  i  v ;  i  ia 
mx'  ,  in  which  ease  trie  cusum  curve  will  fill.  Thus,  in  ; 
wander  up  and  down  about  hr  zero  axis. 

Since  mx’  is  usually  not  known  in  pr.i.-t  i.-e,  .«n  '"itin* 
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Figure  2Qi 


Change  in  construction  operation  as  reflected  in  cumulative  reject 
chart  (schematic). 
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Figure  22  —  Moving  standard  deviation  chart  of  water  content 
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Cumulative  reject  (out-of-specification)  chart  for  compaction 
water  content. 
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Figure  24  —  Distribution  of  compaction  water  content  data  for  tests  1  through 
459. 
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Figure  25  —  Distribution  of  compact 
through  884. 
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PART  VI:  QUALITY  ASSURANCE  BY  ACCEPTANCE  SAMPLING 

The  purpose  of  quality  assurance  is  to  test  fill  as  it  is  placed  and  to 
make  decisions  on  whether  to  accept  or  reject  the  fill  as  conforming  to 
standards.  If  the  fill  is  rejected,  further  compaction  could  be  made,  the  fill 
could  be  removed,  or  some  other  course  of  action  could  be  followed  pursuant  to 
contractual  arrangements  between  contractor  and  owner.  Acceptance  sampling  is 
the  quantitative  tool  used  to  make  the  accept/reject  decision.  The  objective 
of  acceptance  sampling  is  not  to  control  quality,  but  to  make  decisions. 

Structure  of  An  Acceptance  Sampling  Plan 

A  simple  acceptance  sampling  plan  is  structured  in  the  following  way: 

I.  A  random  sample  of  size  n  is  taken  from  the  materials  being 
tested . 

II.  The  results  of  the  n  measurements  (x ■(,..., xn)  are  summarized 

statistically  in  an  index  z.  For  example,  z  might  be  the  sample 
average  (1/n)  Xi . 

III.  The  index  z  is  compared  to  a  critical  value  z*,  and  if  z  lies  on 
the  correct  side  of  z*  the  materials  are  accepted  as 
satisfactory . 

The  questions  in  designing  an  acceptance  sampling  plan  are  how  large  to  make 
the  sample  size  n,  how  to  summarize  the  resulting  data  in  an  index  z,  and  how 
to  select  a  critical  value  z*  such  that  quality  is  assured  without  unduely 
increasing  the  cost  of  construction. 

The  more  stringent  the  acceptance  criteria  become,  the  greater  the 
likelihood  of  rejecting  fill  which  is  in  fact  satisfactory.  The  less 
stringent,  the  greater  the  likelihood  of  accepting  fill  that  is  in  fact  not 
satisfactory.  The  problem  of  acceptance  sampling  is  that,  for  a  given  size 


100 


sample,  reducing  the  likelihood  of  accepting  poor  material  usually  means 
increasing  the  likelihood  of  rejecting  good  materials,  and  vice  versa.  To 
simultaneously  reduce  both  the  likelihood  of  accepting  poor  materials  and  the 
likelihood  of  rejecting  good  materials,  the  sampling  plan  must  be  made  more 
discriminating,  this  usually  increases  inspection  cost. 

Buyer's  Risk  and  Seller's  Risk 

In  specifying  an  acceptance  sampling  scheme  two  risks  are  balanced, 

(a)  The  owner's  (buyer's)  risk  of  accepting  material  of  poor  quality, 
and, 

(b)  The  contractor's  (seller's)  risk  of  rejecting  material  of  good 
quality. 

Decreasing  one  of  these  risks  typically  increases  the  other. 

Test  results  from  an  acceptance  sampling  program  are  variable  whether  the 
fill  is  truly  of  acceptable  quality  or  not.  Because  of  this  variability,  it 
may  be,  for  example,  that  the  lowest  compaction  test  results  on  an  acceptable 
fill  give  lower  dry  densities  than  the  highest  test  results  on  an  unacceptable 
fill. 

The  top  of  Figure  28  shows  a  hypothetical  frequency  distribution  of  test 
results  taken  from  an  acceptable  fill.  Suppose  that  the  criterion  for 
accepting  the  fill  as  meeting  specification  is  that  test  results  be  above  . 
Because  test  results  are  always  variable,  some  fraction  of  the  tests  results 
will  always  fall  below  the  acceptance  criterion  and  thus  lead  to  rejection, 
even  though  the  fill  might  in  fact  be  acceptable.  This  fraction  is 
proportional  to  the  area  under  the  frequency  distribution  to  the  left  of  y^* . 
The  probability  of  the  test  result  lying  beneath  y^* ,  and  therefore  the 


probability  of  improperly  rejecting  an  acceptable  fill,  is  called  the 
seller's  risk. 

The  bottom  of  Figure  28  shows  a  hypothetical  frequency  distribution  of 

test  results  taken  from  an  unacceptable  fill.  Some  fraction  of  these  test 
results  will  always  fall  above  the  acceptance  criterion  yj*  and  thus  lead  to 

the  fill  being  accepted  when  in  fact  it  should  be  rejected.  Ibis  fraction  is 

proportional  to  the  area  under  the  frequency  distribution  to  the  right  of  . 

The  probability  of  the  test  result  lying  above  y^*  and  therefore  leading  to 

acceptance  of  an  unacceptable  fill  is  called  the  buyer1 s  risk. 

For  a  fixed  sampling  plan  there  is  an  explicit  trade  off  between  the 

buyer's  risk  and  the  seller's  risk  in  selecting  the  acceptance  criterion  y^*. 

Higher  values  of  y^*  reduce  the  buyer's  risk  but  raise  the  seller's  risk;  lower 

value  of  y^*  raise  the  buyer's  risk  but  lower  the  seller's  risk.  Ibis  trade 

off  can  be  seen  in  Fig.  28. 

The  buyer's  risk  and  the  seller's  risk  can  be  controlled  simultaneously 
only  by  making  changes  in  the  sampling  plan,  not  just  in  the  acceptance 
criterion.  The  purpose  of  statistical  acceptance  sampling  is  to  allow  the 
buyer's  risk  and  seller's  risk  to  be  quantitatively  determined  for  a  given 
sampling  plan  and  to  be  appropriately  balanced  by  designing  the  sampling  plan. 

Inspecting  for  Fraction  Defective  vs.  Inspecting  for  the  Mean 
Acceptance  sampling  typically  addresses  one  or  both  of  two  aspects  of 
quali ty: 

(a)  The  average  property  of  the  fill,  that  is  the  mean;  or, 

(b)  The  fraction  of  individual  values  within  a  fill  which  are  below  some 
standard,  that  is,  the  fraction  defective. 


Each  aspect  of  quality  may  not  have  the  same  importance  in  a  particular 
application.  For  example,  the  potential  for  internal  erosion  of  a  fill  depends 
on  soil  densities  at  the  least  compacted  places.  Oonversely,  the  strength  of  a 
fill  to  resist  large  slope  instabilities  more  often  depends  on  average  soil 
densities.  Acceptance  sampling  plans  differ  depending  on  which  aspect  of 
quality  is  to  be  assured. 

Operating  Characteristic  Curves 

The  functional  properties  of  an  acceptance  sampling  plan  are  usually 
summarized  by  an  operating  characteristic  or  OC  curve.  The  operating 
characteristic  relates  the  quality  of  the  fill  being  sampled — for  example  its 
mean  density  or  the  fraction  of  the  fill  with  out-of-specification  water 
content — to  the  frequency  with  which  the  sampling  plan  leads  to  a  decision  to 
accept.  As  in  Fig.  29,  the  horizontal  axis  usually  shows  the  actual  fill 
quality,  while  the  vertical  axis  shows  the  probability  of  acceptance.  The 
Buyer's  risk  and  Seller's  risk  are  read  directly  from  the  OC  curve 
corresponding  to  the  definition  of  good  quality  and  poor  quality  materials. 

For  example,  the  probabilities  corresponding  to  the  two  frequency  distributions 
of  Fig.  28  are  shown  as  the  Buyer's  and  Seller's  risk,  respectively,  on  Fig. 

29.  In  principle,  the  better  the  acceptance  sampling  plan,  the  steeper  the  OC 
curve  in  the  vicinity  of  the  contract-specified  quality  of  the  fill.  A  steep 
OC  curve  reduces  both  the  Buyer's  risk  and  the  Seller's  risk. 

The  shape  of  the  OC  curve  depends  on  the  design  of  the  acceptance  sampling 
plan,  and  can  be  used  to  make  economic  decisions  about  the  reasonable  extent 


and  cost  of  sampling.  Usually,  the  easiest  way  to  steepen  the  OC  is  by 


increasing  the  sample  size  and  thus  sampling  cost.  The  remainder  of  Part  VI 


discusses  the  relation  between  a  sampling  plan  and  its  corresponding  operating 
characteristic,  and  how  a  sampling  plan  can  be  designed  to  achieve  a  desired 
OC  curve. 


Acceptance  Sampling  to  Give  Assurance  on  the  Mean 

This  section  considers  acceptance  sampling  plans  the  intent  of  which  is  to 
assure  that  average  properties  of  placed  materials  meet  specification.  TVo 
types  of  speci f ica t ion  are  considered,  single  limits  and  double  limits.  Using 
single  limits  the  concern  is  that  the  average  properties  are,  for  example, 
greater  than  some  specified  value.  For  instance,  average  compacted  dry  density 
is  to  be  greater  than  95”%  standard  or  modified  Proctor  optimum.  Using  double 
limits  the  concern  is  that  the  average  properties  are  between  two  values.  For 
example,  average  compaction  water  content  is  to  be  within  ±2%  standard  or 
modified  Procter  optimum. 

The  sub-sections  first  consider  the  case  of  known  or  specified  material 
variability,  that  is,  known  standard  deviation.  This  case  is  mathematically 
easier  than  the  more  general  case  of  unknown  variability,  and  does  sometimes 
occur  in  practice.  The  more  general  case  of  unknown  variability  is  treated 
afterward. 

■Single  Limit,  Standard  Deviation  Known 

Suppose  specifications  call  for  soil  with  an  averaue  or  mean  compacted  dry 
density  of  mx  =  120  pcf.  Supi»s<»  also  that  the  dry  density  of  the  compacted 
fiLl  is  known  to  have  a  constant  standard  deviation  of  sx'  =  15  pcf.  An 


acceptance  sampling  plan  to  give  assurance  regarding  the  mean  is  constructed 
such  that  material  actually  having  a  mean  of  at  least  120  pcf  (i.e.,  good 
material)  will  be  rejected  no  more  frequently  than  some  fixed  value  a.  As 
before,  a  is  the  Seller's  risk.  Simultaneously  the  sampling  plan  is 
constructed  such  that  material  whose  mean  is  substantially  less  than  120  pcf 
(i.e.,  poor  material)  will  be  accepted  no  more  frequently  than  some  other  fixed 
value  0.  As  before,  8  is  the  Buyer's  risk.  For  the  acceptance  sampling  plan 
to  be  operational,  a  specific  definition  of  what  is  meant  by  "substantially 
less"  must  be  adopted.  In  Fig.  29,  poor  material  is  defined  as  being  an 
average  density  less  than  110  pcf. 

The  procedure  for  acceptance  sampling  with  one  fixed  limit  on  the  mean  is 
the  following 

1 .  Take  a  random  sample  of  n  tests 

2.  From  the  results  x-]  .  .  .xn  calculate  the  mean 

mx--(  1  /n)  Ex-l  . 

3.  Compare  mx  with  a  specified  acceptance  value  m*; 

if  m  >  m*,  then  accept 
if  m  <  m*,  then  reject. 

The  OC  curve  for  a  sampling  plan  regarding  the  mean  shows  the  probability 
of  acceptance  as  a  function  of  the  true  mean  value  of  the  material.  Mx ' ,  as  in 
Fig.  29.  The  OC  curve  is  constructed  by  using  the  standardized  variable  7^, 


The  denomination  in  Bgn.  31  is  the  standard  deviation  of  the  sample  mean 


(1/n)  Zx£  over  repeated  samples  (cf.,  Eqn.  25).  That  is,  the  denomination 
expresses  the  variability  one  naturally  expects  among  different  sets  of  tests. 
The  numerator  is  the  separation  between  the  acceptance  criterion  in  *  and  the 
true  average  quality  of  the  soil  mx'  .  The  variable  Zj^  is  the  number  of 
standard  deviation  separating  m*  for  mx' ,  and  thus  can  be  used  to  calculate  the 
fraction  of  samples  in  which  the  deviation  of  the  sample  mean  for  mx  is  given 
than  m  * *mx  . 

When  the  property  being  tested  has  a  Normal  frequency  distribution,  the 


E 


I 


frequency  distribution  of  zm  over  multiple  samples  is  exactly  Normal.  Yet, 
even  when  the  property  being  tested  does  not  have  a  Normal  frequency 
distribution,  the  frequency  distribution  of  zm  is  still  approximately  Normal. 
The  probability  of  accepting  material  with  actual  mean  mx'  is  found  by 
comparing  zm  with  Table  4  to  find  the  corresponding  frequency  with  which  a 
standard  Normal  variable  exceeds  Zj,, . 

Consider  a  sampling  plan  with  an  acceptance  mean  m*=105  pcf  and  sample 
size  n=6.  Under  this  plan,  n=6  tests  are  made,  the  mean  mx  of  the  results  is 
calculated,  and  if  mx  >  105  pcf  the  material  is  accepted.  If  mx  <  105  pcf  the 
material  is  rejected. 

The  OC  curve  for  this  plan  is  calculated  by  computing  the  quantity  zm  in 
Eqn.  31  and  looking  in  Table  4  to  find  the  probability  of  a  standard  Normal 
variable  having  an  absolute  value  larger  than  Z^  Because  the  standard  Normal 
distribution  is  symmetric  about  ZO ,  the  area  under  the  distribution  above  +Z 
is  the  same  as  the  area  under  the  curve  below  -Z.  Por  example,  if  the  true 
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mean  were  mx  =  120  pcf  and  sx'  =  15  pcf,  then  zm  =  (105  pcf  -120  pcf)/(15 
pcf//6)  =  -2.4.  Thus,  the  probability  of  accepting  good  material  with  the 
specified  mean  density  120  pcf  equals  the  probability  of  a  standard  Normal 
variable  being  algebraically  greater  than  -2.4,  that  is,  about  0.01.  Other 
points  on  an  OC  curve  such  as  that  in  Fig.  29  are  evaluated  by  substituting 
corresponding  values  of  mx'  into  Bgn.  31. 

An  acceptance  sampling  plan  with  regard  to  the  mean  is  designed  by 
specifying  a  Seller's  risk  a  and  a  Buyer's  risk  B.  The  Seller's  risk  is  the 
probability  of  rejecting  fill  which  in  fact  is  of  better  quality  than  some 
decided  upon  acceptable  quality  level  (AQL),  or  "good"  material.  Hie  Buyer's 
risk  is  the  probability  of  accepting  fill  which  is  in  fact  of  poorer  quality 
than  some  decided  upon  unacceptable  quality  level  (UQL),  or  "poor"  material. 

The  AQL  and  UQL  are  engineering  decisions  and  must  be  quantitatively  specified 
to  give  meaning  to  the  notions  of  good  and  poor  quality  material.  Hie  sampling 
plan  is  defined  by  a  sample  size  n  and  an  acceptance  level  m* .  Hie  procedure 
to  find  (n,m*)  is: 

1 .  Specify 

a  =  Seller' s  risk 

S  =  Buye r ' s  risk 

ma  =  Acceptable  quality  level  of  mean  (AQL) 

mu  =  Unacceptable  quality  level  of  mean  (UQL) 

2.  Find  standard  Normal  variables  (Table  4)  with  frequencies  of  not 
being  exceeded  equal  to  ( 1 -a )  and  B, 

z i _a  =  standard  Normal  variable  with  frequency  of  not  being 
exceeded  ( 1 -a  )  . 

zg  =  standard  Normal  variable  with  frequency  of  not  being 
exceeded  B. 


Presumably  for  this  reason.  Normal  distributions  are  common  across  the  broad 
spectrum  of  experimental  science.  In  part  V  of  this  report,  deviations  of 
observed  frequency  distributions  from  Normality  were  used  to  identify  changes 
in  construction  process  and  inspection  procedures. 

Single  Limit,  Standard  Deviation  Unknown 

The  development  of  an  acceptance  sampling  plan  to  assure  the  mean  when  the 
standard  deviation  is  unknown  is  similar  to  the  case  when  the  standard 
deviation  is  known,  except  that  the  index  zm  involving  the  known  standard 
deviation  sx'  is  replaced  by  an  index  t  involving  the  sample  standard  deviation 
sx. 

The  inspection  sampling  procedure  is 

1.  Take  a  random  sample  of  n  tests. 

2.  Calculate  the  mean  and  standard  deviation  of  the  test 
results. 


mx  =  ( 1 /n)  Exi 

-34- 

sx  =  { 1 /n-1 )  E(xi  -  mx)2 

-35- 

3. 

Evaluate  the  sample  statistic 

m  -  m 
^  _  a  x 

-36- 

s  //n 
x 


in  which  ma  =  AQL. 

4.  Fix  an  acceptance  criterion  t*; 

if  t  >  t*  then  accept, 
if  t  <  t*  then  reject. 
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The  OC  curve  for  this  plan  is  a  function  of  both  the  actual  mean  and 


actual  standard  deviation.  The  vertical  axis  of  the  OC  curve  is  the 
probability  of  accepting  the  tested  material.  The  horizontal  axis  is  the 
non-dimensional  quantity  X  =  (ma ' -mx ' )/sx ' ,  which  involves  the  specified  AQL 
and  both  the  real  mean  and  the  real  standard  deviation  cf.,  plot  in  S.W.. 

The  OC  curve  is  calculated  in  a  manner  similar  to  that  when  the  standard 
deviation  is  known,  but  whereas  the  variation  of  zm  across  different  samples 
can  be  approximated  by  a  Normal  frequency  distribution,  the  variation  of  t 
across  different  samples  --  at  least  for  small  n  (Say  n  >  20)  —  is  wider  than 
for  zm  and  must  be  approximated  by  a  so-called  student  -t  frequency 
distribution,  The  variation  in  t  is  wider  than  in  zm  because  the  sample 
standard  deviation  varies  somewhat  from  sample  to  sample.  As  n  gets  larger, 
the  variation  of  sx  about  sx'  becomes  smaller,  and  the  student  t  distribution 
approaches  a  Normal  distribution. 

To  design  an  acceptance  sampling  plan  for  the  mean  when  the  standard 
deviation  is  unknown  the  procedure  is: 

1.  Specify, 

a  =  Seller' s  risk 

B  =  Buyer's  risk 

ma  =  Acceptable  (mean)  Quality  Level  (AQL) 

mu  =  Unacceptable  (mean)  Quality  Level  (UQL) 

2.  Make  a  rough  estimate  of  sx' 


3. 


Compute  X  = 


m  -  m 
a _ u  . 

s 


-38- 


This  is  the  value  of  A  when  the  actual  mean  mx'  equals  the  UQL  mu' 


4.  Make  a  rough  estimate  of  n  from  Pig.  32. 

Material  having  the  value  of  A  from  step  3  should  be  accepted 
only  g  fraction  of  the  time.  Find  n  in  Fig.  32  providing  B 
probability  of  accepting  material  of  quality  A 

5.  Find  the  acceptance  criterion  t*  corresponding  to  a  frequency  of 
not  being  exceeded  (1-a)  from  Table  6,  using  v  =  n-1  degrees  of 
freedom . 

6.  Specify  sampling  plan  by 

n  =  sample  size 

t*  =  acceptance  criterion 


t  =  - 77 -  =  test  statistic.  -39- 

s  //n 
x 


Plate  4  shows  the  design  of  a  sampling  plan  for  the  same  condition  as  in 
Plate  3,  but  that  the  standard  deviation  is  not  known.  The  effect  of  not 
knowing  the  standard  deviation  in  this  case  is  that  the  sample  size  must  be 
increased  by  one  test,  from  9  to  10,  to  obtain  the  same  precision  in  the  OC 
curve . 


Double  Specification  Limits,  Standard  Deviation  Known 

Certain  material  properties,  as  for  example  compaction  water  content,  have 
specification  limits  both  above  and  belcw  their  target  value.  Soil  moisture 
should  be  within  some  ±  interval  of  optimum,  say,  no  wetter  than  +2%  of 
standard  or  modified  optimum  Proctor  and  no  dryer  than  -2%.  *in  acceptance 
sampling  plan  with  double  specification  limits  intends  to  assure  that  a 
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material  property  is  within  the  defined  interval. 

An  acceptance  sampling  plan  with  double  limits  is  designed  by  specifying 
two  acceptance  bounds.  If  the  sample  average  lies  between  these  bounds,  the 
lift  is  accepted.  If  the  sample  average  lies  outside,  the  lift  is  rejected. 

The  bounds  are  chosen  to  conform  to  specific  values  of  the  Buyer's  risk  and  the 
Seller's  risk. 

If  the  variability  of  the  soil  properties  as  measured  by  the  standard 
deviation  is  known,  then  the  variability  of  the  sample  average  of  n  tests  from 
one  sample  to  another  is  also  known  (i.e.,  sm=sx'//n).  As  before,  if  the  soil 
properties  are  assumed  to  have  a  Normal  frequency  distribution,  the  variation 
of  the  sample  average  also  has  a  Normal  distribution.  Even  if  the  soil 
properties  are  not  Normally  distributed,  the  distribution  of  the  sample  average 
is  usually  still  approximately  Normal. 

Let  the  target  value  or  acceptable  quality  level  of  the  average  soil 
properties  be  ma.  If  indeed  the  average  soil  property  is  ma,  the  sample 
average  of  n  test  results  will  vary  about  ma  as  shown  in  Fig.  33.  This 
sample-to-sample  variability  of  mx  is  centered  on  ma  and  has  standard  deviation 
sx//n. 

Let  my*  and  my*  be  the  upper  and  lower  acceptance  limits  on  the  sample 
mean  mx.  If  mx  is  greater  than  my*  or  less  than  my*  the  lift  is  rejected.  The 
Seller's  risk  a  is  the  frequency  with  which  the  sample  mean  mx  lies  outside  my* 
and  my*  when  in  fact  the  true  mean  is  ma .  That  is,  the  Seller's  risk  is  the 
shaded  areas  in  under  the  frequency  distribution  of  mx  in  Fig.  33.  Each  tail 
area  has  frequency  (i.e.,  probability)  a/2. 


Lee  UQL(j  and  UQLl  be  the  upper  and  lower  unacceptable  quality  levels.  If 
the  actual  average  soil  property  lies  just  outside  the  UQLp  or  UQLL,  there  is 
still  a  chance  that  sample  variability  will  allow  the  measured  sample  mean  to 
lie  inside  the  range  (mL*,  my*),  and  thus  lead  to  the  lift  being  improperly 
accepted.  This  frequency  is  the  Buyer's  risk  8»  The  sampling  variability  of 
mx  for  two  lifts  which  have  true  means  equal  to  UQLp  and  UQLl  are  shown  in  Fig. 
34.  The  frequency  (i.e.,  probability)  with  which  the  sample  mean  from  these 
soils  lies  within  ..he  interval  (mL*,  my*)  is  shown  by  the  shaded  areas  under 
the  respective  frequency  distributions.  Each  tail  area  equals  8,  the  Buyer's 
risk. 

To  design  an  acceptance  sampling  plan  on  the  mean  with  double 

specification  limits,  two  constraints  must  be  satisfied,  the  Seller's  risk  and 

the  Buyer's  risk.  TVo  parameters  can  be  adjusted,  the  sample  size  n  and  the 

location  of  the  acceptance  limits  mp*  and  m^/  •  The  sample  size  controls 

the  width  of  the  frequency  distribution  of  mx,  in  that  the  standard  deviation 
of  mx  equals  sx'//n;  while  the  limits  mg*  and  mL*  control  where  the  frequency 

distributions  are  cutoff  to  yield  a  and  3. 

From  Table  4,  the  tail  area  under  a  Normal  frequency  distribution  can  be 
related  to  numbers  of  standard  deviation  on  either  side  of  the  mean.  Let  Zp  be 
the  number  of  standard  deviations  below  which  the  area  under  the  Normal 
frequency  distribution  is  (1-p)  (i.e.,  Zp  is  the  standard  Normal  variable  which 
has  probability  p  of  not  being  exceeded).  For  example,  from  Table  4,  zo.97b  = 

+  1.96,  and  Zq.q2S  ~  -1.96.  Then,  Figs.  33  and  34  lead  to  four  relationships 
from  which  an  acceptance  sampling  plan  can  be  designed: 


As  an  example,  consider  an  inspection  plan  for  compaction  water  content  in 
which  the  Seller's  risk  and  Buyer's  risk  were  set  at  a=0  .05  and  B=O»10, 
respectively.  The  target  value  of  average  water  content  is  Proctor  optimum, 
and  intolerable  deviation  from  the  target  has  been  decided  to  be  +3%  water 
content.  Assume  that  from  project  records  the  standard  deviation  were  known  to 
be  about  1.5%.  For  these  conditions,  Bqns.  40  to  43  become 


r  *  n  wif  iru 


mu  "  0 
1  .5//n 


* 

m  -  O 

Lj 

1 .5//n 

mu  “  3% 
1.5//n 


m  -(-3%) 

Lj 

1 .  5//n 


=  -1.96 


=  -1.282 


1.282 


Solving  the  first  two  equations  simultaneously  gives. 


*  * 

mU  =  ~  mL 


Solving  the  first  and  third  equation  simultaneously  gives, 


n  =  [(-!_)  (1.5)  (1.96  +  1 . 282) ] 2 


=  2.62, 


or  rounding  off,  n=3.  Putting  n= 3  into  the  equation  for  Seller's  risk  gives 
mu*  =  1.7%,  mL  =  -1.7%.  Putting  n=3  into  the  equation  for  Buyers  risk  gives 
mu*  =  1.9%  and  mL  =  -1.9%.  Choosing  ±1.8%  as  the  acceptance  limits  gives  a 
Seller's  risk  of  a=0 . 38  (i.e.,  less  than  5%)  and  a  Buyer's  risk  of  8=0.08 


(i.e.,  less  than  10%).  Ihe  OC  curve  for  this  plan  is  shown  in  Fig.  35. 

Another  example  is  given  in  Plate  5. 

Double  Specification  Limits,  Standard  Deviation  Unknown 

When  the  standard  deviation  is  unknown  the  procedure  for  specifying  an 
acceptance  sampling  plan  is  much  the  same  as  when  the  standard  deviation  is 
known,  except  that  the  sample  standard  deviation  sx  replaces  the  known  standard 
deviation  sx'  in  Bqns.  40  to  43,  and  the  Student-t  distribution  (Table  6) 
replaces  the  Normal  distribution  (Table  4). 

As  for  the  case  of  a  single  specification  limit,  the  test  statistic  is, 


t 


m  -  m 


x _ a 

s  //n 
x 


-50- 


For  sample  sizes  above  about  n=20  these  modifications  are  unnecessary  because 
the  sample  standard  deviation  sx  is  sufficiently  close  to  the  actual  standard 
deviation  sx' . 

The  inspection  sampling  procedure  is, 

1 .  Take  a  sample  of  size  n 

2.  Calculate  the  mean  and  standard  deviation  of  the  test 
results, 

mx  =  (1/n)£  xi  -51- 

sx  =  ( 1 /n-1 )  £(x i  -  mx)2  -52- 

3.  Evaluate  the  sample  statistic 
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Fix  an  acceptance  criterion  t* ; 

if  t  <  t*  then  accept 
if  t  >  t*  then  reject. 


As  in  the  case  of  a  single  specification  limit,  an  inspection  sampling  plan 
with  double  limits  is  designed  by  specifying  a  sample  size  n  and  a  criterion 
t* .  An  initial  guess  at  sx'  is  made,  and  Fig.  32  is  used  to  estimate  a  sample 
size  n  based  on  the  quantity. 


|UQL  -  m 


in  which  ma  is  the  target  soil  property  and  UQL  is  either  the  upper  or  lower 
unacceptable  quality  level.  This  assunes  that  UQL|j  and  UQLl  are  symmetrically 
placed  about  the  target  ma.  See  Duncan  (1974)  or  Grant  and  Leavenworth  (1972) 
for  asymmetric  cases.  Unacceptable  materials  at  either  the  UQLy  or  UQLl  should 
be  accepted  only  with  frequency  B*  Thus,  knowing  the  Buyer's  risk  6  and  the 
number  of  standard  deviations  \  separating  the  UQL  from  ma,  an  initial  sample 
size  can  be  chosen  from  Fig.  32.  Using  UQL=±3%,  ma=0%,  and  sx  =1.5%,  as 
before,  Fig.  32  leads  to  n  ~  4. 

The  acceptance  criterion  t*  is  found  from  a  table  of  the  Student' s-t 
frequency  distribution  (Table  6).  This  table  provides  the  frequencies  with 
which  given  values  of  the  test  statistic  of  Eqn.  50  are  exceeded  due  to  random 
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sampling  variations  when  in  fact  the  soils  being  inspected  are  of  target 


quality  ma.  Because  both  unacceptably  high  and  unacceptably  low  values  will  be 
rejected,  the  Seller's  risk  is  the  sum  of  the  frequencies  with  which  the  test 
index  of  Bqn.  50  lies  above  +  t*  and  below  -t*  .  Thus,  t*  is  set  so  that  the 
tail  areas  on  either  side  each  have  probability  a/2.  For  the  Student' s-t 
frequency  distribution  these  tail  areas  depend  on  the  sample  size  taken,  though 
the  so-called  degrees  of  freedom  v=n-1 •  For  these  conditions.  Table  6  gives  a 
t*  value  of  3.25. 


5 .  Specify  Sampling  Plan: 


a)  Take  sample  of  4  tests 

b)  Calculate  sample  mean  mx  and  sample  standard  deviation  sx 

c)  Calculate  the  test  index 


m  -  0 
x _ 

t  s  //4 
x 

d)  If  -3.18  <  t  <;  +3.18,  then  accept  lift. 
If  t  <  -3.18  or  t  >  +3.18,  then  reject. 


-55- 


Acceptance  Sampling  for  Fraction  Defective 
The  following  section  considers  the  case  in  which  an  inspection  sampling 
plan  is  employed  to  assure  that  the  fraction  of  out-of-specification  material 
in  a  compacted  fill  is  within  tolerable  limits.  Such  plans  are  generally 
called  acceptance  sampling  for  fraction  defective. 

If  the  soils  data  display  a  Normal  or  bell-shaped  frequence  distribution, 
there  is  an  exact  mathematical  relationship  between  the  mean  and 
deviation  of  the  data  on  the  one  hand,  and  the  fraction  defective  on  the  other. 
This  is  shown  schematically  in  Fig.  36. 


Setting  the  lower  limit  of  acceptable  values  or  the  UQL  at 


L  =  lower  limit  of  acceptability 


a  standardized  deviate  z^  is  defined  as  the  number  of  standard  deviations  sx' 
separating  the  mean  of  the  data  mx  from  the  lower  limit  L, 


m '  -  L 

x 


-5  6- 


s ' 

x 


If  the  data  are  normally  distributed,  Zl  is  uniquely  related  to  the  fraction 
defective,  as  shown  in  Fig,  36.  Numerical  values  of  this  relationship  are 
found  in  Table  4  or  can  be  approximated  by  Eqn.  24.  Fig.  36  illustrates  that 
the  higher  the  mean  and  the  lower  the  standard  deviation,  the  lower  the 
fraction  defective.  An  acceptance  sampling  plan  for  fraction  defective  is 
structured  in  the  following  way: 


1.  Test  a  random  sample  of  size  n  to  obtain  the  data  x-|,..,,xn. 

2.  From  the  results,  calculate  a  sample  mean  mx,  sample  standard 
deviation  sx,  and  test  index 


m  -  L 
x 


-5  7- 


Depending  on  the  specific  problem,  the  formula  for  z  may  vary 
somewha  t. 


Compare  the  computed  value  of  z  with  a  critical  value  z* : 
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f.  v. 

k.-.V.  f't. 


/, 


A 


If  z  >  z* ,  then  accept. 

If  z  <  z* ,  then  reject. 

The  choice  of  n  and  z*  defines  the  performance  properties  of  the  sampling  plan. 
These  parameters  are  usually  chosen  to  satisfy  specified  levels  of  Buyer's  risk 
and  Seller's  risk. 

Operating  Characteristic  for  Fraction  Defective  Sampling 

The  operating  characteristic  or  OC  curve  summarizes  the  discriminatory 
power  of  an  acceptance  sampling  plan.  The  OC  curve  shows  how  the  probability 
of  accepting  a  lift  or  other  quantity  of  material  varies  as  a  function  of  the 
quality  of  the  material  being  inspected.  For  plans  aimed  at  fraction  defective 
the  OC  curve  relates  probability  of  acceptance  to  the  fraction  defective  in  the 
lift. 

Consider  an  acceptance  sampling  plan  for  percent  compaction  specified  by: 

n  =5 

L  =  95%  maximum  Proctor  density 

z*  =  1.645  (i.e.,  5%  of  the  soil  less  than  95%  max.  density). 
Presume  the  standard  deviation  is  known  to  be  sx’  =  2%.  Under  this  plan  5 
tests  are  made.  The  average  of  the  tests  mx  is  compared  to  L  through  z  = 

(mx  -  L)/sx’.  If  z  >  1.645  the  material  is  accepted;  if  z  <  1.645  the  material 
is  rejected.  Ihe  OC  curve  for  this  plan  relates  the  probability  of  accepting 
the  material  to  the  actual  fraction  of  the  lift  compacted  to  less  than  95% 
Procter  maximum. 

For  Normally  distributed  material  with  known  standard  deviation  there  is  a 
unique  relation  between  the  fraction  defective  and  the  mean.  For  sx'  =  2%  and 
L,  =  95%  Table  4  is  used  to  find  the  following  relations: 


Fraction  defective 

P' 

0  .05 

0.10 

0.15 

0.20 

0.25 

Mean 

mx' 

98.  3 

97.6 

97. 1 

96.7 

96.4 

Thus,  the  horizontal  axis  of  the  OC  curve  can  be  expressed  either  as  actual 
fraction  defective  or  as  actual  mean. 

For  a  given  fraction  defective  or  given  mean,  the  probability  of  accepting 
the  material  equals  the  probability  that  the  test  result  z  is  greater  than  z*  = 
1.645.  This  probability  can  be  determined  by  noting  that  z  is  itself  Normally 
distributed.  With  L  and  sx'  fixed,  z  depends  only  on  the  mean  mx  of  the  test 
results.  When  sampling  from  a  Normally  distributed  population,  the  frequency 
distribution  of  the  sample  mean  is  also  Normal  (Part  II).  Thus,  the 
probability  that  z>1.645  is  found  by  calculating  the  mean  and  standard 
deviation  of  z  and  referring  to  Table  4. 

The  mean  and  standard  deviation  of  z  are  found  by  the  method  described  in 


Part  II, 


Table  4  is  entered  by  calculating  the  number  of  standard  deviation  of  z 
separating  mz  from  the  acceptance  criterion  z*=1.645.  The  corresponding  number 
on  the  vertical  axis  is  the  probability  of  rejection  (i.e.,  the  tail  area  of 
the  Normal  curve,  or  Pr{z<z*}).  The  probability  of  accepting  the  material  is 
the  complement  of  this  number, 

Pr {accepting^  =  i  _  Pr| rejecting! .  -60- 

This  procedure  is  illustrated  in  Plate  6. 

The  entire  OC  curve  is  found  by  calculating  the  probability  of  accepting 
the  material  for  various  values  of  actual  fraction  defective.  For  the  sampling 
plan  above,  the  full  OC  curve  is  shown  in  Fig.  37.  If  none  of  the  material  is 
defective  the  probability  of  accepting  is  1.0,  and  as  the  actual  fraction 
defective  increases  (i.e.,  as  the  mean  of  the  material  decreases)  the 
probability  of  accepting  goes  down. 

Single  Limit  with  Known  Standard  Deviation 

The  main  question  in  designing  an  acceptance  sampling  plan  is  to  decide 
upon  a  sample  size  n  and  an  acceptance  criterion  z* .  These  choices  dictate  how 
the  plan  performs  with  respect  to  Buyer's  risk  and  Seller’s  risk.  Let  the 
probability  of  improperly  accepting  unsatisfactory  material,  the  Buyer's  risk 
be, 


Buyer's  risk  =  8;  -61- 

let  the  probability  of  improperly  rejecting  satisfactory  material,  the 


Seller' s  risk  be. 


Seller's  risk  =  a. 


-62- 


An  OC  curve  is  defined  by  specifying  two  points  through  which  it  passes. 
For  this  purpose  the  Buyer  must  specify  a  maximum  fraction  defective  that  he 
considers  tolerable  and  which  would  be  accepted  under  the  plan  only  some 
fraction  (?  of  the  time.  This  poor  quality  material  as  a  fraction  defective  is 
denoted  Py' .  At  the  same  time,  a  target  or  desired  quality  level  is  specified 
which  would  be  accepted  at  least  (1-a)  fraction  of  the  time.  This  good  (i.e., 
acceptable)  quality  as  a  fraction  defective  is  denoted  pa' .  Hie  OC  curve  can 
be  made  to  pass  through  the  two  points  (pa',1-a)  and  (Pu',8)  by  adjusting  the 
sample  size  n  and  acceptance  criterion  z* . 

For  example,  consider  that  acceptable  material  has  pa'=0.01  fraction 
defective  and  an  unacceptable  material  has  py'sO.IO  fraction  defective.  To  fix 
the  two  points  of  the  OC  curve  specified  by  the  Buyer's  risk  and  the  Seller's 
risk,  the  first  task  is  to  calculate  the  corresponding  averages  ma'  and  mu' 
which  would  give  fractions  defective  of  pa'=C.01  and  pu'=0.10,  respectively. 
From  Table  4,  the  area  under  the  Normal  curve  below  -2.33  standard  deviations 
from  the  mean  equals  0.01,  and  the  area  below  -1.28  standard  deviation  equals 
0.10,  Tius,  an  acceptable  soil  having  pa'=O,01  and  standard  deviation  sx'=2% 
would  ha/e  a  mean, 

ma  =  L  +  2.33  sx' 

=  95%  +  2. 33(2%)  -63- 

=  99.7%; 
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and  an  unacceptable  soil  having  py'^O.IO  would  have  a  mean, 
mu  =  L  +  0 . 84  sx' 

=  95%  +  1.28(2%)  -64- 

=  97.6%. 

The  test  index  z  is  calculated  from  Bqn.  56.  IXie  to  random  sampling 
variability,  the  value  of  z  varies  from  one  sample  of  n  tests  to  another  even 
for  the  same  soil.  This  sampling  variability  can  be  characterized  by  a  mean  mz 
and  standard  deviation  sz  for  each  of  the  soils  above.  Specifically,  for  the 
acceptable  quality  soil. 


m 

z 


99.7-95 

2 


2.  35 


-65- 


sz  =  1//n 


-66- 


For  the  unacceptable  quality  soil,  mz  =  1.28  and  sz  =  1//n.  These  are  the 
means  and  standard  deviations  that  the  test  statistic  z  would  have  if  the 
actual  soil  being  tested  were  just  at  the  edge  of  being  acceptable  or  just  at 
the  edge  of  being  unacceptable,  respectively. 

The  Buyer's  risk  and  Seller's  risk  specify  target  probabilities  of 
accepting  the  two  types  of  soil  above  when  using  the  acceptance  sampling  plan. 
For  acceptable  soil  Pr{z<z*}  =  a;  for  unacceptable  soil  Pr{z>z*}=8.  This  gives 
two  equations.  Again  from  Table  4,  for  Pr{z<z*}  =  a  =  0.05  the  mean  of  z  must 


be  1.645  standard  deviation  larger  than  z* , 


2.33  -  1 . 645 // n 


-67- 


=  z* 


For  Pr{z>z*}  =  6  =  0.10,  the  mean  of  z  must  be  1.28  standard  deviations  smaller 
than  z* , 


mzu  "*■  1  •  28  Sj;u  -  z* 

1.28  +  1 . 28//n  =  z*  . 

Bqns.  67  and  68  are  solved  simultaneously  to  give, 


-68- 


n  =  7.79  +  say,  8  -69- 

z*  =  1.75  -70- 


An  example  of  the  acceptance  sampling  plan  is  specified  as  shown  in  Plate  7. 

The  design  of  an  acceptance  sampling  scheme  may  be  accomplished  more 
quickly  by  algebraically  solving  for  the  sample  size  n  and  acceptance  criterion 
z*.  Define, 


z  1  -a 

z  1  -6 


standardized  Normal  variable  for  which  the 
probability  of  not  being  exceeded  (Table  6)  is  1-a. 

standardized  Normal  variable  for  which  the 
probability  of  not  being  exceeded  (Table  6)  is  1-  6« 

za  =  standardized  Normal  variable  for  which  the 

probability  of  not  being  exceeded  (Table  6)  is  1-pa' 

zu  =  standardized  Normal  variable  for  which  the 

probability  of  not  being  exceeded  (Table  6)  is  1-Pu' 

For  a  single  criterion  acceptance  sampling  plan  having  parameters  (pa',a)  and 


(pu',8)  the  sample  size  and  acceptance  criterion  are. 


To  summarize,  the  procedure  for  designing  a  single  limit  acceptance 


sampling  plan  is: 

1.  Select  a  Seller's  risk  a,  and  a  Buyer's  risk  B* 

2.  Select  acceptable  quality  level  pa'  and  unacceptable  quality  level 
Pu'  • 

3.  Find  values  for  standard  Normal  variables  corresponding  to  1 -a ,  1- 
(1-pa'),  and  (  1  — pu '  )  probabilities  of  not  being  exceeded  (z-|_a,  Zf 
za '  zu ) • 

4.  Calculate  the  sample  size  by 


n 


z  -t-  z 
1-a  1-6  1 2 

z  -  z  ‘ 
a  u 


-73- 


5.  Calculate  the  acceptance  criterion  by 


* 

z 


z  z  +  z  z 

a  1-8  u  1-a 

z  *-  z 

1-a  1-B 


-74- 


6 


Plot  the  (X  curve 


Single  Limit  with  Unknown  Standard  Deviation 


Usually  the  standard  deviation  of  the  material  property  being  tested  is 
unknown.  The  only  information  about  the  standard  deviation  comes  from  the  data 
themselves  in  the  form  of  the  sample  standard  deviation, 


/ 


=  /  — T  (x.  -  m  ) ' 
n-1  L  1  x 


-75- 


If  the  sample  size  is  large  (n>20),  the  sample  standard  deviation  will  be  close 
to  the  real  standard  deviation  and  the  assumption  of  known  standard  deviation 
can  be  made  with  neglible  error.  If  the  sample  size  is  not  large,  a  slight 
modification  to  the  foregoing  procedure  must  be  made. 

When  the  standard  deviation  is  unknown  the  quantity  z  is  calculated  using 
the  sample  standard  deviation  sx. 


z  = 


m  -  L 
x 


-76- 


Whereas,  when  the  standard  deviation  is  known  the  quantity  z  has  a  Normal 
distribution,  when  the  sample  standard  deviation  is  substituted  for  the  real 
standard  deviation  the  calculated  value  of  z  has  more  variability.  Now  the 
denominator  as  well  as  the  numerator  will  vary  from  one  sample  to  another.  The 
frequency  distribution  of  z  takes  on  the  slightly  broader  shape  of  the 


Student  -t  distribution. 


The  procedure  for  designing  sampling  plans  and  calculating  OC  curves  when 
the  standard  deviation  is  unknown  is  the  same  as  when  the  standard  deviation  is 
known,  with  the  exception  that  tables  of  the  Student  -t  frequency  distribution 
rather  than  the  Normal  distribution  are  used.  Areas  under  the  Student  -t 
distribution  for  the  standardized  case  of  zero-mean  and  unit  standard  deviation 
are  given  in  Table  6.  Note,  unlike  the  Normal  distribution,  Student  -t 
depends  on  the  sample  size  n.  As  n  becomes  large  the  shape  of  the  Student  -t 
approaches  the  Normal  distribution. 

Convenient  approximations  for  sample  size  and  acceptance  criterion  when 
the  standard  deviation  is  unknown  are  (Wallis,  1947), 


z  z  +  z  z 
*  a  1-B  u  1-a 

z  = 


Z.  +  z_  „ 
1-a  1-B 


-77- 


n 


(l  +  z*2/2) 


(V -*■■!—  )2 


z  -  z 
a  u 


-78- 


Thus,  when  the  standard  deviation  is  not  known  a  larger  sample  must  be  taken  to 
get  the  same  OC  curve.  The  sample  size  must  be  larger  by  the  factor  (1+z*2/2). 

The  example  of  the  previous  section  is  recalculated  in  Plate  8,  now 
relaxing  the  assumption  that  the  standard  deviation  is  known.  Ttie  OC  curve  can 
be  calculated  approximately  but  acceptably  by  assuming  z  to  be  Normally 
distributed  with  a  standard  deviation  equal  to  sx  [l/n  +  (  z*2 /2n )  }  ^ /2  .  The 
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approximate  OC  curve  is  shown  in  Fig.  39.  Thus,  to  calculate  the  real  fraction 


defective  p'  corresponding  to  a  given  probability  of  acceptance  q  (i.e,  to  plot 


the  OC  curve) ,  first  the  corresponding  standardized  Normal  deviation  Zq  is 


taken  from  Table  4.  Next,  Zg  is  increased  by  the  factor  ( 1 /n  +  ( z*2 /2n)  ) 1 /2  . 


Then  a  corresponding  zD'  is  calculated  as 


Zpi  =  z*  -  zq  (l/n  +  (z*2/2n))1/2 


Then  Table  4  is  used  to  determine  p' .  For  example,  in  Plate  7,  z*  =  1.63  and 


n  =  6. 


Double  Specification  Limits 


The  preceding  plan  pertains  to  the  case  of  one  specification  limit.  Fbr 


example,  dry  density  should  be  at  least  95%  standard  or  modified  Procter 


maximum.  When  deviations  in  either  direction  are  important  the  plan  must  be 


modified.  For  a  lower  limit  of  acceptability  L  and  an  upper  limit  U,  the 


minimum  fraction  defective  occurs  when  mx  lies  halfway  between  L  and  U;  that 


is,  when  the  limits  are  symmetric  about  the  mean.  In  this  case, 


m  -  L 
x 


U  -  m 


and 


B 


/  V  >  .*  V 


The  fraction  defective  equals  the  area  under  the  Normal  curve  outside  ±z,  or 


twice  the  fraction  defective  read  from  Table  4.  Note  that  z  depends  only  on 
the  upper  and  lower  limits  IJ,  L  and  on  sx '  .  It  does  not  depend  on  mx'.  Hius, 
if  sx '  is  known,  the  first  step  is  to  assume  that  the  acceptable  fraction 
defective  is  yreater  than  the  area  under  the  Normal  curve  outside  +  z.  This  may 
be  done  without  samplinq,  and  indicates  whether  the  variability  of  the 
construction  process  reflected  m  sx '  is  so  larqe  as  to  preclude  any 
possibility  of  the  tested  soil  beinq  found  acceptable. 

Presuininy  that  z  from  Dpi.  Ho  is  sufficiently  larqe  that  rejection  is  not 
inescapable,  the  fraction  defective  will  dejiend  on  both  Zq  and  zg,  and  the 
acceptance  criterion  must  ne  based  on  both.  In  concept,  this  is  done  by 
summing  the  fraction  defective  beneath  L  and  the  fraction  defective  above  U  and 
comparinq  that  sum  to  the  criterion  M.  However,  a  simpler  procedure  can  be 
developed  by  considerinq  the  operatinq  character istic  curve  of  the  sampling 
plan. 


Double  Limits,  Standard  Deviation  Known 

The  fraction  defective  for  double  specification  limits  is  that  proportion 
of  the  area  under  the  frequency  distribution  of  the  material  property  which 
lies  either  below  a  lower  specification  limit,  ^ 

L  =  lower  specification  limit. 


or  above  an  upper  specification  limit, 

U  =  upper  specification  limit. 

For  constant  standard  deviation  the  fraction  defective  is  minimized  when  the 
mean  mx'  lies  halfway  between  L  and  U.  In  this  case, 

-(L-mx')/sx'  =  (U-mx')sx'  =  <U-L)/2sx'  -83- 

Thus,  a  quick  check  should  be  made  to  see  whether  a  material  can  possibly  meet 
the  fraction  defective  double  specification  standard  by  finding  the  area  under 
the  Normal  frequency  distribution  outside  ±  z  =  (U-L)/2sx'.  If  this  area  is 
greater  than  the  acceptable  fraction  defective  pa'  no  sampling  plan  alone  will 
assure  quality.  The  construction  process  must  be  changed  to  make  the  material 
more  uniform  and  thus  reduce  sx' . 

In  the  general  case  for  double  specification  limits,  an  acceptance 
sampling  plan  to  assure  fraction  defective  follows  the  following  procedure: 

1 .  Take  a  sample  of  size  n 

2.  From  the  results,  calculate  the  sample  mean 
mx  =  ( 1 /u) Exi 

3.  Compute  the  quantities 


4.  Specify  an  acceptance  criterion  z*  :  If  zL  >  z*  an^  Zjj  >  z*  , 
then  accept  otherwise,  reject. 

The  problem  with  double  specification  limits  is  determining  z* .  In  the 

case  of  single  specification  limit  z*  was  determined  from  areas  under  the 
Normal  frequency  distribution  to  one  side  of  a  specification  limit.  In  the 

double  specification  case  z*  must  be  determined  from  the  sum  of  the  areas  above 

U  and  below  L. 

Consider  the  problem  of  acceptance  sampling  for  compaction  water  content. 
The  target  value  is  Proctor  optimum  water  content.  The  upper  speci f ication 
limit  is  U  =  +2%  optimum;  the  lower  specification  limit  is  L  =  -2%.  Presuming 
the  standard  deviation  of  water  content  to  be  1%,  the  limits  are 


U-L  +2  -  (-2) 

f  z  =  +  -  -  +■ 


2s 


2(1  ) 


+  2 


-86- 


That  is  L  and  U  are  4  (i.e.,  + 2 )  standard  deviations  apart.  From  Table  4  the 
area  under  the  Normal  curve  beyond  z=2  is  0.02.  Hius,  the  lowest  possible 
fraction  defective  would  be  twice  0.02  or  about  4%.  The  fraction  defective  for 
values  of  the  mean  other  than  that  halfway  between  U  and  L  are  shown  in  Table 

8. 

Presume  f;>r  sake  of  example  that  the  acceptable  quality  level  or  AQL 
expressed  as  a  fraction  defective  were  pa '  =  0.10.  Tbat  is,  the  lift  would  be 
considered  acceptable  Lf  at  least  got  of  the  soil  had  a  compaction  water 
content  between  I,  -  -2%  Proctor  optimum  and  U  -  +2%.  From  Table  8  (by 


interpolating  my 


values)  any  lift  with  an  average  water  content  bet  :en  -0.7% 


and  +0.7%  would  be  acceptable,  for  the  fraction  of  any  of  these  lifts  with 


water  contents  outside  ±2%  would  be  less  than  0.10.  At  mx'  =  -0.7%  the 
fraction  below  -2%  is  0.096  and  the  fraction  above  +2%  is  0.004.  The  sum  is 
0.10.  Similarly  but  in  reverse  at  mx'  =  +0.7%,  the  fraction  above  +2%  is  0.096 
and  the  fraction  below  -2%  is  0.004.  The  double  limit  specification  can  thus 
be  met  by  combining  two  single  limit  tests  designed  such  that  the  acceptable 
fraction  defective  in  each  is  reduced  from  pa'  =0.10  to  pa'  =  0.096.  One 
applies  on  the  upper  limit  side,  the  other  applies  on  the  lower  limit  side. 

The  design  for  these  two  plans  is  exactly  as  discussed  before,  and  is  carried 
out  in  Flate  9. 

Double  Limits,  Standard  Deviation  Unknown 

The  problem  of  designing  an  acceptance  sampling  plan  for  fraction 
defective  with  double  sped  f  ication  limits  and  unknown  standard  deviation  is 
less  easily  solved  than  the  single  limit  problem.  In  particular,  with  double 
limits  the  shape  of  the  OC  curve  depends  on  how  the  fraction  defective  is  split 
between  the  upper  and  lower  tail  of  the  distribution.  However,  the 
availability  of  statistical  tables  and  graphs  designed  expressly  for  the 
purpose  (US  DOD  Military  Standard  4 1 J )  jreatly  simplifies  the  task.  For  the 
purpose  of  acceptance  sampling  of  engineered  fills,  the  jraphs  of  Fig.  41  and 
Fig.  42  provide  sufficient  accuracy. 

The  procedure  begins  is  for  the  single  limit,  unknown  standard  deviation 
case.  Buyer’s  and  Seller's  risk  •»  and  are  specified,  and  acceptable  and 
unacceptable  fractions  defective  pa '  and  p,,'.  ttpis.  77  and  "M  arc  used  to 


estimate  a  sample  size  n  and  an  acceptance  criterion  z*. 


From  the  estimates  of 


n  and  z*  the  quantity 

1  - 

Y  =  - 


* 

z  /n 
n  -  1 

2 


-87- 


is  calculated  and  used  to  enter  the  absicca  of  Fig.  41.  On  the  ordinate,  and 
corresponding  to  the  appropriate  value  of  n,  an  allowable  fraction  defective  M 
is  read. 

The  test  procedure  is  implemented  by  taking  a  sample  size  n,  calculating 
the  sample  mean  mx  and  sample  standard  deviation  sx,  and  then  computing  the 
test  indices 


z 


L 


-88- 


z 


U 


-89- 


From  Fig.  42,  for  the  appropriate  value  of  n,  estimated  fractions  defective 
corresponding  to  zL  and  zg  are  read  of  as  pg  and  pg,  respectively.  These  are 
summed  to  obtain  an  estimate  of  the  total  fraction  defective,  which  is  in  turn 
compared  to  M  to  decide  whether  to  accept  or  reject  the  lift.  If  pg  +  pg  <  M, 
then  the  lift  is  accepted;  otherwise,  the  lift  is  rejected.  An  example  is 
given  in  Plate  10.  The  OC  curve  for  this  procedure  is  approximately  the  same 


as  in  the  single  limit  case,  using  the  same  values  of  a,0,  Pa'*  an<3  Pu' 


Table  8  —  Fraction  defective  for  double  sampling  limits. 


mx 

L-m  v 
— -$• 
sx 

U~fx 

sx 

PU 

P=PL+PU 

+  2.0 

4.0 

0.0 

0.50 

0.50 

+  1  .5 

3.5 

0.5 

- 

0.31 

0.31 

+  1  .0 

3.0 

1  .0 

0.001 

0.16 

0.17 

+0.5 

2.5 

1 .5 

0.006 

0.07 

0.0  86 

0  .0 

2.0 

2.0 

0.02 

0.02 

0.04 

-0.5 

1  .5 

2.5 

0.07 

0.006 

0  .086 

-1  .0 

1  .0 

3.0 

0.16 

0  .001 

0.17 

-1  .5 

0.5 

3.5 

0. 31 

- 

0. 31 

-2.0 

0 

4 

0.50 

- 

0.50 

Note:  range  U-L  is  kept  constant 


SUBJECT:  Acceptance  sampling  plan  to  assure  mean  value  of 

compacted  dry  density;  standard  deviation  known. 


SUBJECT:  Acceptance  sampling  plan  to  assure  mean  value  of  compacted  dry 

density;  standard  deviation  unknown. 


I.  PROBLEM:  Design  acceptance  sampling  plan  to  assure  mean  value  of  compacted 

dry  density. 

II.  SOLUTION: 

1 .  Specify, 

a  =0.05 

6  =0.10 

ma  =  AQL  =  1  20  pcf 
mu  =  UQL  =  110  pcf 

2.  Estimate  Standard  Deviation 

sx'  a  10  pcf. 

3.  Estimate  Sample  Size 


ma  mu  1 20  -  1 10  pcf . 

A  =  7  =  To~^f  •  1,0 

s 

x 

For  A=1 .0  the  probability  of  accepting  should  be  @=0.10.  Fig.  32 
shows  that  n=10  is  the  approximate  sample  size  providing  this 
probability  of  accepting. 

4.  Findt* , 


From  Table  6  the  value  of  t  which  is  exceeded  (1-a)=  0.95  fraction 
of  the  time  is  t*=-1.83. 
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PLATE  5 


SHEET  1/2 


SUBJECT:  Design  acceptance  sampling  plan  to  assure  average  compaction  water 

content  within  double  speci f ica tion  limits;  standard  deviation  is 
known . 


I.  PROBLEM: 

Design  acceptance  sampling  plan  to  assure  average  compaction 

water  content  is  within  ±2%  Procter  optimum.  Standard  deviation  is  known 

to  be  sx  =2%. 

II.  SOLUTION: 


1 .  Specify: 

a=0  .05 
6=0  .05 
5x' 


Target  water  content  =  Procter  optimum. 
UQLy  =  +2%  Procter  optimum. 

UQLl  =  -2%  Procter  optimum. 


2.  Write  eguations  for  Seller's  risk  and  Buyer's  risk: 


rau 

-  m 

a 

* 

S 

X 

//n 

m* 

Li 

-  m 

a 

i 

S 

X 

//n 

m* 

-  u2Lu 

s 

X 

//n 

m* 

Li 

"  u^l 

*1  -a/2 


"a/2 


mu  - 


2//n 


2//n 


m*  -  2% 


=  1.96 


=  -  1.96 


'1-8 


2//n 


m*  -  (-2%) 


1 .645 


=  1.645 


s  //n 
x 


2// n 


I  1  ) 


f-4 


■ . 1 

■>  • 


>1 


0 

r  si 


% 


8 

■’.N 

* 
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SUBJECT: 


Design  acceptance  sampling  plan  to  assure  average  compaction  water 
content  within  double  specification  limits;  standard  deviation  is 
known . 


Solve  equa tions  simultaneously: 


a) 

from 

Bqns.  1 

and  2, 

mU* 

=  -mL* 

b) 

from 

Bqns.  1 

and  3 , 

n  » 

1  3 

c) 

from 

Bqns.  1 

,2,  3,4, 

mtJ* 

=  -mL* 

=  1 

.09 

Spe< 

=ify 

sampling 

plan : 

a) 

Take 

sample 

of  size  n 

=  1  3 

wa  ter 

contents 

b) 

Calculate  mean  of  sam 

pie 

mx  = 

( 1 /n) 

c) 

If 

-1  .09% 

<  mx  <  +1 

.09% 

Procter  optimum,  then 

accept  lift. 

d) 

If 

^  1  • 

09%  or  mx 

<  1 

.09% 

Procter  optimum. 

then  reject 

lift. 

NOTE:  The  relatively  large  sample  size  and  tight  acceptance  limits  for  this 

inspection  plan  are  caused  by  the  large  variability  of  the  soil  relative  to  the 
unacceptable  quality  limits  of  +2%  Proctor  optimum  water  content. 


SUBJECT:  Operating  characteristic  (OC)  curve  for  fraction  defective  sampling 

typical  calculation;  standard  deviation  known. 
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PLATE  7 


SHEET  1/3 


SUBJECT:  Acceptance  sampling  plan  for  fraction  defective;  standard  deviation 

known . 


PROBLEM: 

Design  an  acceptance  sampling  plan  for  fraction  defective 
%-compaction,  when  the  standard  deviation  is  known. 

SOLUTION: 

1.  Specify:  Buyer's  Risk  a  =  0.05 
Seller's  Risk  8  =  0.10 

Acceptable  fraction  defective  pa '  =  0.01 
Unacceptable  fraction  defective  py'  =0.10 
Standard  deviation  sx '  =2%  Procter  optimum 
Lower  Limit  of  acceptable  compaction  L=95%  Procter 
optimum 

2.  Calculate  average  percent  compaction  corresponding  to  pa '  a 


Acceptable  compaction: 


99.7% 


Unacceptable  compaction; 


Define  test  index: 


=  97.6% 


SUBJECT:  Acceptance  sampling  plan  for  fraction  defective;  standard  deviation 


PLATE  8 


SHEET  1/2 


SUBJECT:  Acceptance  sampling  plan  for  fraction  defective;  standard  deviation 

unknown . 


PROBLEM:  Design  an  acceptance  sampling  plan  for  soil  compaction  with  the 

properties 

a  =0 .05  pa’  =0.01 

8  =0.10  pu'  =  0.20 

L  =  95%  optimum  Procter 

sx  =  Unknown. 

SOLUTION: 

1 .  Find  standard  Normal  variables  corresponding  to  ct ,  8  /  Pa  1  >  Pu 1  • 


From  Table  4:  a  ♦  z1-a  =  1*645 

8  2 1  =  1.28 

pa'  +  za  =  2.33 

Pu'  -*■  zu  =  0.84 

2 .  Calculate  sample  size  n  arid  acceptance  criterion  2*. 


z  zM  „  +  z  z„ 

*  a  1-0  u  1-a 
z  = 


z  +  z 
1-a  1-8 


=  (2.33)  (1.28)  +  (  0.84)  (  1.645) 

1.28  +  1.645 


=  1  .5 


n 


(1 


(1 


+ 


1.492 

2 


z  +  z 
,  1-a  1-  8}2 

z  -  z 
a  u 


1 .645  +  1.28 
2.33  -0.84 


8.13  +  8 


PLATE  8 
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SUBJECT: 


Acceptance  sampling  plan  for  fraction  defective;  standard  deviation 
unknown . 


V’ 
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SUBJECT:  Acceptance  sampling  plan  to  assure  fraction  defective  with  double 

specification  limits;  standard  deviation  known. 


PROBLEM: 

Design  an  acceptance  sampling  plan  to  assure  fraction  defective  on 
the  basis  of  compacted  water  content.  Assume  SD  known  to  be  1%  optimum 
Proctor . 


SOLUTION: 


Specify, 


U  =  +2%  Procter  optimum 
L  =  -2%  Procter  optimum 


2.  Determine  values  of  mx'  such  that  the  sum  of  fraction  defective  above 


U  and  fraction  defective  below  L  equals  the  AQL,  Pa'=0.10 


from  TabLe  8  (by  interpolation)  : 
Jan  PU  PL  P  TOTAL 


PLATE  9 


SHEET  2/2 


Acceptance  sampling  plan  to  assure  fraction  defective  with  double 
specification  limits;  standard  deviation  known. 


Evaluate  n  and  z* , 


From  Eqn.  71  n 


/Z 1  -a+  Z 1  -8  -j? 

^  z  -  z  ' 
a  u 

+  1 .65+1 . 28  s2 


-  (fr 


30  -0  . 5  3 


+  14.5  +  15 


Prom  Rgn.  72  z*  = 


Za  2  1  -8  +  Zu  2 1 -a 
Z1-a  +  Z1-B 


=  (1  .  30)  (1  .  28)  +  (0 .53)  (1 .65  )  =  0.87 

( 1 . 65 )  +  ( 1 . 28) 


Specify  acceptance  sampling  plan 

a)  Take  random  sample  of  size  15 

b)  Calculate  sample  mean  mx,  and 

U  -  mx 

zu  =  - ~ 


test  indices, 


z 


L 


m  -  L 
x 


i 


s 

X 


c)  Compare  zu  and  zL  with  z*  =0.87, 
If  zu  and  zL  >  z* ,  then  accept 
If  zu  and  zL  <  z* ,  then  reject 

PC  Curve  is  shown  in  Fig.  40. 


PLATE  10 


SHEET  1/2 


SUBJECT:  Acceptance  sampling  plan  to  assure  fraction  defective  with  double 

specification  limits;  standard  deviation  unknown 


PROBLEM : 

Design  an  acceptance  sampling  plan  to  assure  fraction  defective  on 
the  basis  of  compaction  water  content.  Standard  deviation  is  unknown 

SOLUTION: 

1 .  Specify, 


a  =  0.10 

8  =  0.10 

Pa'  =  °-05 

pu'  =  0.30 


U  =  +2%  Procter  optimum 
L  =  -2%  Procter  optimum 


2.  Determine  standard  Normal  deviates  corresponding  to  a, 8, Pa  and 


a  =0.10 
8  =0.10 
pd'  =  0.05 
Pu'  =  0 . 30 


Zi_a  =  +1.28 
=  +1.28 
za  =1.65 
z,,  =+0.5  3 


Use  Bqns •  77  and  78  to  estimate  n  and  z*, 


*9,,.  1  -a  +  Z 1  -8  1 2  1.02?  ,+1.28+1.28^2 

+z*2/2  l-Y—; - -)  =  M+— >  (usV-O-Vr^ 

a  u 

Za  Z1-8  *  Zu  ?'l-u  _  (  1 .65  )  ( 1  .  28)  +  (0  .53)  (1  .  28)  = 

z  +  z,  „  "  (  1  .  28)  +  { 1 . 28) 

1 -a  1-8 


7.94  +  8 


=  1.09 


Estimate  allowable  fraction  defective  M 


Compute  quantity 


1  -  z *  /n 


(n-1  )  ^  1-( 1 . 0  9) /8 | (8-1  )  _ 

2  2 


Read  M  fr'>n  Fig.  41  for  y=0.28,  n  =  8  ♦  M  =  0.15 


Acceptance  sampling  plan  to  assure  fraction  defective  with  double 
specification  limits;  standard  deviation  unknown 


Specify  test  procedure 


a.  Take  a  random  sample  of  size  n=8 

b.  Compute  sample  mean  mx,  and  standard  deviation  sx,  calculate 
test  indices, 


c.  From  Fig.  42  read  values  of  pa  and  pu  corresponding  to  Zjj  an 
ZL 

d.  Compare  p^  +  py  to  M, 

If  PL  +  Pu  <  M*  then  accept. 

If  Pl  +  Pu  <  M»  then  reject. 

The  PC  curve  for  this  plan  is  shown  in  Fig.  43 


Figure  33  —  Frequency  distribution  of  the  sample  average  mx  for  a  sample  of 
size  n. 


Probab 


•ji  '  *  WJ  ’  •  V  V  v.v 


Probability  of  Acceptance 


Figure  40  —  Operating  characteristic  curve  for  sampling  plan  aerivec  m  Plate  g 


—  Chart  for  determining  fraction  defective  from  z  (from  Duncan, 


For  Standard  Deviation  Plans  Take  abscissa  =  - ^  — — 


For  Average  Range  Plans  lake  abscissa  =  - — 


Figure  42 


Chart  for  determining  maximum  allowable  fraction  defective 
(from  Duncan,  1974). 
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